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ABSTRACT 


This  Final  Technical  Report  on  Contract  DAAH01-7J  -C-1192  discusses 
investigations  of  tracking  errors  associated  with  a  particular  monopulse 
phasea-array  radar  system,  called  the  Experimental  Array  Radar  (EAR).  The 
free-space  performance  of  the  EAR  is  first  analyzed  and  degradations  of 
performance  due  to  particular  target  Doppler  shifts  and  multipath  returns 
are  noted.  A  computer  analysis  of  the  tracking  performance  of  the  EAR  is  then 
presented,  and  it  is  noted  that  there  are  significant  tracking  errors  due  to 
multipath  returns  for  targets  at  altitudes  less  than  1200  feet.  An  analysis 
of  the  ef fectiveness  of  frequency-agile  operation  in  reducing  the  effects 
of  these  multipath  returns  is  presented,  which  shows  that  the  available 
frequency-agility  bandwidth  of  a  modified  EAR  would  significantly  reduce 
multipath-induced  tracking  errors,  particularly  for  higher  antenna  and  target 
locations.  Improved  Mil  filters  which  provide  more  nearly  uniform  responses 
for  a  wide  range  of  Doppler  frequencies  are  then  discussed  and  a  number  of 
representative  improved  filter  responses'  presented. 
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I .  INTRODUCTION 


This  Final  Technical  Report  describes  v;ork  performed  for  the  U.  S.  Army 
Missile  Command  under  Contract  DAAH01-71-C-1192  from  May  1971  through  May  1972. 
The  investigations  carried  out  during  this  period  involved  analyzing  the 
Experimental  Array  Radar  (EAR)  being  constructed  by  the  Missile  Systems 
Laboratory,  and  developing  methods  for  improving  its  performance. 

These  investigations  were  initially  based  on  wgrk  performed  for  the 
Missile  Command  under  Contract  DAAH01-70-C-0535  [l]  ,  under  which  Georgia 
Tech  examined  the  effects  of  polarization  agility  on  monopulse  radar  angle 
tracking.  The  techniques  and  target  models  developed  under  this  earlier 
Missile  Command  contract  formed  the  basis  for  a  substantial  portion  of  the 
work  performed  under  the  current  contract. 

The  initial  phase  of  the  program  was  concerned  with  preliminary  analyses 
of  the  performance  of  the  EAR  system  and  identification  of  areas  of  marginal 
performance.  A  modified  C-Iland  radar  at  Georgia  Tech  was  used  to  estimate 
anticipated  radar  cross-section  of  targets  for  incorporation  into  this 
preliminary  performance  analysis.  Results  of  this  preliminary  analysis  are 
presented  in  Section  II,  and  indicate  deterioration  in  the  performance  of  the 
EAR  when  the  target  has  a  relatively  low  Doppler  frequency  and/or  when  strong 
multipath  returns  are  being  received  by  the  system. 

The  effect  of  strong  multipath  returns  on  the  tracking  performance  of 
the  EAR  when  tracking  low-flying  targets  is  analyzed  in  some  detail  in  Section 
III.  The  applicability  of  frequency  agility  to  the  EAR  system  for  the  purpose 
of  minimizing  these  multipath-induced  errors  is  discussed  in  Section  IV.  The 
analyses  in  both  of  these  chapters  leans  heavily  on  a  computer  prediction  of 
the  EAR  performance  which  included  realistic  representation  of  antenna 
beamshapes,  null  positions,  and  sidelobe  levels. 

The  design  of  optimum  MTI  processors  for  use  in  the  EAR  system  is  treated 
in  some  detail  in  Section  V.  Previous  work  in  digital  processing  is  reviewed 
and  improved  procedures  for  design  of  processors  having  optimum  responses 
for  a  range  of  Doppler  shifts  are  presented. 

Conclusions  and  recommendations  resulting  from  this  study  are  presented 
in  Section  VI. 


Numbers  in  brackets  refer  to  References  in  Section  VII. 
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II.  PRELIMINARY  INVESTIGATION  OF  EAR  PERFORMANCE 


The  initial  step  in  the  research  program  was  a  preliminary  system 
analysis  to  define  the  performance  of  the  EAR  and  to  focus  attention  to  those 
areas  where  performance  is  marginal.  First, the  free-space  signal-to-noise 
ratio  was  calculated  for  the  EAR  system,  and  the  effect  of  thermal  noise  on 
angle-tracking  accuracy  was  analyzed.  Next,  the  effect  of  multipath  returns 
entering  the  antenna  sidelobes  (the  high  angle  tracking  case)  was  considered, 
and  finally  signal-to-clutter  ratios  were  calculated  for  the  EAR  system 

A.  Free-Space  Signal-to-Noise  Ratio 

The  system  parameters  assumed  in  the  preliminary  analysis  of  the  free- 
space  signal-to-noise  ratio  for  the  EAR  system  were: 

P  =  peak  transmitted  power  =  75  kW; 

G  =  antenna  gain  (assumed  to  be  same  for  both  transmit  and 
receive  modes)  =  25  dB; 

_2 

X  =  wavelength  =  5.45  x  10  meters  (5.5  GHz);  and 


N  =  noise  power  =  -93  dBm  (10-dB  noise  figure,  10-MHz  bandwidth). 

The  free-space,  single-pulse,  signal-to-noise  ratio  (S/N),  is  given  by 


S/N  = 


2  2 

PGlo 

t _ 

3  4 

N(4tt)  R 


2.245  x  10  a 
R4 


where 

R  =  range  in  meters,  and 

a  =  target  radar  cross-section  in  square  meters. 

The  average  signal-to-noise  ratio  (averaged  over  all  Doppler  frequencies 
and  relative  phases)  is  unaffected  by  MTI  processing.  This  can  be  seen 
readily  by  considering  that  the  average  power  gain  of  a  two-channel  MTI 
processor  is  six,  and  that  six  noise  samples  are  added  during  processing,  so 
that  the  output  signal-to-noise  is  unaffected.  In  a  single-channel  processor, 
the  average  gain  is  reduced  to  three  (due  to  the  necessity  of  averaging  over 
all  relative  phases)  but  now  only  three  noise  pulses  are  added,  so  the  output 
signal-to-noise  ratio  remains  unchanged. 

When  processing  in  triplets,  as  in  the  EAR  system  [2*1,  the  48  pulses  on 
a  target  produce  16  independent  samples.  Approximately  10  dB  should  be  added 
to  the  signal-to-noise  ratio  calculated  earlier  to  account  for  effects  of 
integration  of  these  16  samples.  Figure  1  shows  the  average  free-space  signal- 
to-noise  ratio  as  a  function  of  range  for  several  values  of  target  idar 
cross-section. 
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Average  signal-to-noise  ratio  for  the  EAR  system  as  a 
function  of  range.  These  values  should  be  modified  by 
the  response  at  the  specific  target  Doppler  frequency 
obtained  from  Figure  2. 
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The  average  signal-to-noise  ratio  is  given  by  the  above  expression; 
however,  the  response  at  a  particular  Doppler  frequency  may  be  considerably 
different  from  this  average  value.  The  frequency  response  for  the  conventional 
three-pulse  MTI  filter  is  derived  in  Appendix  A  and  plotted  in  Figure  2. 

Figure  2  has  been  normalized  so  that  the  signal-to-noise  ratio  for  a  given 
target  may  be  determined  by  adding  the  response  at  the  desired  Doppler  frequency 
from  Figure  2  to  the  values  read  from  Figure  1.  For  very  slow  targets  and  for 
targets  near  the  blind  speed,  substantial  reductions  in  performance  are 
possible.  The  design  of  improved  MTI  filters  having  more  nearly  uniform 
response  is  discussed  in  Section  V. 


B.  Tracking  Errors  Due  to  Glint  and  Thermal  Noise 

The  rms  tracking  error  due  to  thermal  noise  has  been  derived  by  a 
number  of  authors  to  be  [3] 


k  v/Bt  (S/N) 
m 


where 

at  =  rms  angle-tracking  error, 

0  =  3  dB  antenna  beamwidth, 

k  =  difference-channel  error  slope, 
m 

B  =  i-f  bandwidth, 

T  =  pulse  length,  and 
S/N  =  signal-to-noise  power  ratic  , 

If  one  assumes  Bt  =  2,  k  =1.57,  and  9=2°,  which  might  be  reasonable  choices, 
then 

0.9 

at-  =  -  > 

c  /s7FT 

where  at  is  expressed  in  degrees. 

This  rms  tracking  error  in  degrees  is  plotted  as  a  function  of  range 
in  Figure  3  for  a  number  of  target  radar  cross-sections.  The  10  dB  of  inte¬ 
gration  gain  included  in  Figure  1  is  also  incorporated  into  the  results  shown 
in  Figure  3.  The  signal-to-noise  ratio  used  was  the  average  over  all 
Doppler  frequencie.;  at  some  particular  frequencies,  the  performance  may  be 
appreciably  worse  due  to  the  frequency  response  of  the  OTI  processor. 

The  treatment  of  rms  error  due  to  target  glint  may  be  considerably 
simplified  by  characterizing  the  target  by  an  ''..ms  effective  length,"  Leff. 
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Figure  2.  Response  of  a  conventional  three-pulse  canceller  having 
filter  tap  weights  of  1,  -2,  and  1.  The  response  has  be 
normalized  so  the  signal  power  averaged  over  all  Doppler 
frequencies  is  0  dB. 


RMS  Tracking  Error  (Degrees) 


1  2  3  4  5  6  7  8  9  10  15  20 

Range  (km) 


Figure  3.  Average  tracking  errors  due  to  thermal  noise  and  glint  as 

Functions  of  range  for  the  EAR  system.  Integration  improvement 
in  signal-to-noise  ratio  of  10  dB  is  included. 
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The  rms  angular  error  in  degrees  as  a  function  of  range  is  also  plotted  in 
Figure  3  for  targets  having  several  values  of  L  Results  of  analyses 
performed  during  earlier  contracts  at  Georgia  Teen  and  surveys  of  the 
available  literature  indicate  that  an  effective  length  of  from  5  to  10  meters 
in  azimuth  is  representative  of  values  to  be  expected  from  many  targets  of 
interest.  The  effective  length  in  elevation  will  probably  be  somewhat  less; 
reasonable  values  probably  lie  between  1  and  3  meters. 

Figure  3  indicates  that  the  region  of  maximum  average  angle  tracking 
accuracy  lies  roughly  between  3  and  7  km,  depending  upon  the  specific  choices 
of  target  radar  cross-section  and  effective  target  length  assumed. 

C.  Radar  Angle  Tracking  Errors  Due  to  Multipath  Returns 

The  preceding  analyses  indicate  the  magnitude  of  errors  to  be  expected 
when  targets  are  being  tracked  under  free-space  conditions.  Additional 
tracking  errors  are  introduced  by  reflections  from  the  surface  of  the  earth. 

Two  different  methods  of  analysis  have  been  used  to  study  these  tracking 
errors.  Barton  [3]  discusses  the  case  where  the  interfering  signal  is  considered 
to  be  noiselike  with  random  phase,  and  with  amplitude  determined  both  by  the 
sidelobe  ratio  for  the  direction  of  arrival  of  the  multipath  signal  and  by  the 
reflection  coefficient  of  the  reflecting  surface.  Another  method  of  analysis 
assumes  the  multipath  signal  to  be  a  return  from  a  second  target  whose 
amplitude  is  given  by  the  product  of  the  sidelobe  ratio  and  the  reflection 
coefficient  of  the  reflecting  surface.  The  errors  are  then  those  of  this 
two-scatterer  target  as  their  relative  phase  changes.  Comparing  results  of 
these  two  analyses,  plotted  in  Figure  4,  shows  them  to  be  in  substantial 
agreement. 

For  high  elevation  angles,  an  average  sidelobe  level  of  25-dB  and  a 
surface  reflectivity  of  0.5  results  in  a  multipath  signal  approximately  31  dB 
down.  From  Figure  4,  this  corresponds  to  an  error  of  approximately  0.018 
beamwidths  or  0.04  degrees;  such  errors  will  only  limit  the  tracking  accuracy 
for  targets  having  large  cross-sections  near  the  region  of  maximum  accuracy. 

While  this  limits  the  maximum  accuracy  attainable  for  1-10  m  cross-section 
targets  for  ranges  approximately  5-10  km,  the  problem  becomes  worse  for 
low-angle  targets,  wnere  the  image  target  is  more  strongly  illuminated  and 
the  multipath  return  represents  a  target  which  is  stronger  and  less  random 
in  nature,  A  more  detailed  computer  analysis  of  radar  tracking  performance 
in  this  region  is  discussed  in  Section  III. 

D.  Clutter  Effects  on  the  EAR 

The  limitations  on  clutter  cancellation  with  digital  MTI  as  implemented 
in  the  EAR  system  appear  to  be  due  to  quantization  errors  in  the  analog-to- 
digital  conversion  process,  rather  than  the  actual  clutter  residue  associated 
with  the  frequency  response  of  the  filter  and  the  width  of  the  clutter  spectra. 
The  calculations  in  Appendices  A  and  B  indicate  the  clutter  residue  produced 
by  a  conventional  three -pulse  canceller  processing  return  from  wooded  hills  in 
a  20-knot  wind  is  approximately  -83  dB,  while  the  limitation  on  cancellation 
due  to  analog-to-digital  quantization  error  for  an  eight-bit-plus-sign 
conversion  is  approximately  -53  dB. 
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RMS  Error  (Bearawidchs) 


The  received  signal-to-clutter  power  ratio  (S/C)  before  MTI  processing 
may  be  approximated  by 


S/C 


ct  S  S 
rt  rg 

o 

a 


x 


X 


_2 

TC 


where 


cr  =  target  radar  cross-section, 

S  =  gain  in  the  sum  channel  on  reception  for  target, 

S  =  gain  in  the  sum  channel  on  reception  for  reflected  energy, 

o°  =  radar  cross-section  per  unit  area  of  clutter, 

R  =  range, 

6.  =  azimuth  3-dB  beamwidth, 

A 

t  =  pulse  length,  and 
c  -  velocity  of  light. 

For  the  vertical  array,  where  the  broad  azimuth  beamwidth  (approximately  30°) 
results  in  large  clutter  power  and  the  broad  transmitting  pattern  is  assumed 
to  illuminate  the  clutter  and  the  target  with  approximatley  equal  intensity, 
the  S/C  ratio  is  a  function  of  elevation  angle.  However,  if  elevation 
greater  than  approximately  2  ,  the  gain  for  the  clutter  signal  will  be 
approximately  the  average  gain  in  the  sidelobes  of  the  vertical  array.  For 
simplicity,  assume  that  this  gain  is  approximately  -25  dB.  Values  of  S/C 
ratio  may  now  be  calculated  and  then  modified  by  the  clutter  cancellation  of 
the  system. 

A  plot  of  the  average  signal-to-clutter  ratio  after  cancellation  for  a 
one-microsecond  pulse  and  various  rations  of  a/a  valid  for  elevation  angles 
greater  than  approximately  2  is  presented  in  Figure  5.  These  values  must 
be  modified  by  the  response  of  the  MTI  processor  for  the  particular  target 
velocity. 

2  o 

To  illustrate  the  use  of  Figure  5,  consider  a  =  1  m  and  a  -  -20  dB  or 
a/a°  =  100  tn  =  20  dBsm.  The  average  signal-to-clutter  ratio  is  approximately 
43  dB  at  a  range  of  10  km.  However,  the  frequency  response  of  the  MTI 
processor  must  be  added  to  this  value;  since  the  processor  response  may  be 
substantially  smaller  than  -20  dB,  the  value  obtained  from  Figure  5  may  be 
modified  substantially.  This  reduction  of  signol-to-clutter  ratio  is  another 
manifestation  of  nonuniform  Doppler  responses  produced  by  conventional  MTI 
processors,  and  illustrates  the  need  for  improved  MTI  filter  responses  as 
discussed  in  Section  V. 
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III.  MULTIPATH  EFFECTS  ON  EAR  ANGLE  TRACKING  ACCURACY 


A.  Introduction 

In  the  previous  section,  it  was  noted  that  multipath  reflections  returned 
to  the  antenna  through  the  low-amplitude  sidelobes  of  the  EAR  antenna  do  not 
limit  system  performance  except  for  the  case  of  large  targets  near  the  region 
of  maximum  accuracy.  However,  when  target  elevation  is  reduced,  multipath 
returns  become  stronger,  errors  are  less  random  in  nature,  and  tracking 
accuracy  lessens;  consequently,  multipath  signals  in  this  situation  can 
seriously  limit  EAR  system  performance. 

Characteristics  of  the  target,  reflecting  surface,  tracking  antenna,  and 
radar  data  processor  all  influence  the  nature  and  severity  of  multipath- 
induced  tracking  errors.  A  computer  program  was  written  to  assist  in 
determining  the  effects  of  these  factors  on  tracking  performance  of  the  EAR 
in  multipath  situations.  The  program  permits  calculation  of  elevation-angle 
tracking  error,  and  signal  amplitude  and  phase  variations,  and  allows  imple¬ 
mentation  of  various  processing  techniques. 

The  geometry  of  the  target  trajectory,  calculation  of  received  signals 
and  indicated  errors,  antenna  pattern  generation,  program  organization,  and 
results  of  this  computer  analysis  are  discussed  in  the  following  paragraphs. 

B.  Geometry  of  the  Target  Trajectory 

The  analysis  assumes  a  target  flying  at  constant  altitude  on  a  radial 
path  toward  the  antenna  as  shown  in  Figure  6.  Target  altitude  is  denoted  by 
H  and  the  antenna  is  located  at  height  H&.  The  surface  is  assumed  to  be 
fiat,  so  that  simple  specular  reflection  occurs,  and  to  hav<  a  reflection 
coefficient  p  (in  most  of  the  analyses  p  was  assumed  to  be  real).  The  angles 
necessary  to  define  the  geometry  are  indicated  in  Figure  1  as  follows: 

THT  Angle  between  antenna  and  target  measured  from  horizontal 

THIM  Angle  between  antenna  and  target  image  measured  from 

horizontal 

THANT  Pointing  angle  of  antenna  normal  measured  from  horizontal 

THSCAN  Angle  at  which  antenna  beam  is  scanned  from  antenna  normal 

THBEAM  Pointing  angle  of  antenna  beam  measured  from  horizontal 

TH1  Angle  between  antenna  beam  and  target  measured  from  antenna 

base 

TH2  Angle  between  antenna  beam  and  target  image  measured  from 

antenna  beam 


Preceding  page  blank 
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Note:  All  angles  are  measured 
positive  in  the  CCW 
direction.  Dotted  arcs 
denote  angles  which  are 
negative  for  the  con¬ 
figuration  shown. 


Figure  6.  Antenna/target  geometry. 
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V.  < 


C.  Calculation  of  Received  Signals  and  Indicated  Errors 

In  the  EAR  system,  the  signals  are  transmitted  through  the  horizontal 
array  and  received  through  both  the  horizontal  and  vertical  arrays;  therefore, 
different  antenna  gains  apply  on  transmission  and  reception  when  tracking  with 
the  vertical  array.  With  reference  to  Figure  6,  let  S  and  D  represent 
the  sum-pattern  gain  and  difference-pattern  gain,  respectively^  of  the  phased- 
array  antenna  upon  reception  in  the  a  direction;  similarly  S  will  denote 
the  transmitting  antenna  gain  in  the  a  direction.  The  subscript  a  will  denote 
either  the  direction  of  the  direct  ray  to  the  target  (  a  -  t)  or  the  ray 
reflected  from  the  ground  to  the  target  (a  =  g) .  When  ref lectio rs  are  present, 
the  total  signal  received  by  the  antenna  will  be  a  combination  o.t  signals 
which  traverse  four  paths.  The  first  path  is  that  directly  to  the  target 
and  back  to  the  antenna,  the  second  and  third  are  paths  having  one  reflection 
from  the  ground  surface,  and  the  fourth  involves  two  reflections  from  thx' 
ground  surface.  Thus,  the  total  signal  received  at  the  antenna  (for  the  t-um 
channel)  is  given  by 


S  =  S  e-JkR  S  „  e-jkR  +  S  e-jkE  p  S  e 
tt  rt  tt  ^  rg 


-jk(R  +  6R)  s  e-jkR 


rt 


-jk(R*6R)  e‘j(kR  +  6R) 

tg  rg 

where  p  is  the  ground  reflection  coefficient  and  6R  is  the  range  difference 
between  direct  and  reflected  paths  to  the  target.  The  decrease  in  signal 
level  due  to  the  increase  in  range  by  an  amount  6R  is  neglected  as  a  second- 
order  effe^^  Simplifying  and  changing  the  phase  reference  to  delete  the 


factor  e 


there  results: 


S  =  S.  .S  (SS  +  S  S  )  pe' jk6R  +  p2S.  S  e"j2k6R  . 
tt  rt  tt  rg  tg  rt  v  tg  rg 

In  a  similar  fashion,  the  difference  signal  strength  is  found  to  be: 

D  ■  ht\  +  »sCtY'lk6R  +  pStSV'3kSR  +  ^psY'3"611  • 

The  indicated  error,  9  ,  may  then  be  calculated  from  the  ratio  of  the 
difference  signal  to  the  sum  signal  as 

8e  =  K  Im  (-§-)  , 

where  K  is  a  constant  chosen  to  give  unity  error  slope  at  boresight  and  is 
a  function  of  antenna  pattern  beamwidth  and  beamshape.  A  typical  plot  of 
indicated  angular  error  versus  actual  angular  error  is  presented  in  Figure  7 
for  K  =  0.86.  Because  the  error  functions  become  less  sensitive  as  the  antenna 
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Actual  Angular  Error  (degrees) 


is  scanned  (due  to  increasing  antenna  beamwidth)  it  may  be  advantageous  to 
r;;ake  K  a  function  of  scan  angle  so  the  curves  for  all  scan  angles  have 
approximately  unity  slope  at  the  origin.  Curves  for  this  case  are  shown  in 
Figure  8  where  K  =  0.845  +  0.335  (0  scan/30). 

D.  Antenna  Pattern  Generation 

For  the  system  simulation,  it  was  decided  to  use  a  simple  scheme  that 
would  give  an  adequate  representation  of  the  antenna  patterns  to  be  approxi 
mated.  The  characteristics  of  the  antenna  pattern  used  for  the  computer 
analysis  of  vertical  array  are  shown  in  Table  I,  and  compared  with  measured 
values  given  in  parentheses: 

Table  I  -  Antenna  Pattern  Characteristics 

Parameter 
Sum  3  dB  BW 

Gain  Reduction 

Sidelobe  Level  -27  dB  (-28)  -30  dB  (-28) 

Peak  Difference  -2  dB  (-2.7)  -2.5  dB  (-2.7) 

Pattern  Level 
(referred  to  peak  of 
sum  patterns) 

The  patterns  were  generated  using  assumed  aperture  illuminations  of 
the  f^rm 

f(x)  =  A  cos  , 

which  result  in  far-field  patterns  having  the  form 

/  \  _  I  r.  cos  u 
g(u)  -  k.F  2  > 

TT  2 

2  •  u 

where 

u  =  k2?  sin  9^ 

To  generate  the  sum  and  difference  patterns,  0.  is  assigned  values 
9  and  9  +  0  ,  where  9  is  the  offset  angle  of  the  beam.  The  two  patterns 
thus  created  are  adde8  and  subtracted  to  produce  sum  and  difference  far-field 
patterns  (see  Figure  9).  In  order  to  cause  the  patterns  to  vary  with  scan 
angle,  both  0q  and  F  are  made  to  be  functions  of  scan  angle,  the  specific 
function  being  determined  empirically  to  produce  the  desired  pattern  variation. 

The  transmitted  signal  is  radiated  through  the  horizontal  array  which 
has  very  nearly  constant  gain  over  a  wide  range  of  elevation  angles.  The 
elevation  pattern  for  this  antenna  was  simulated  as  being  that  of  an  isotropic 
radiator. 
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E.  Program  Organization 

The  computer  program  used  for  analysis  (Appendix  E)  consists  of  a  large 
main  program  which  performs  input  and  output  tasks  with  several  subroutines 
to  do  specialized  computations.  All  parameters  are  initialized  in  the  main 
program;  this  program  controls  the  execution  of  the  various  subroutines. 

The  subroutine  ERRO,  for  example,  provides  for  calculation  of  the  indicated 
angular  error  between  antenna  and  target  and  returns  this  information  to  the 
main  program  for  antenna  beam-pointing  calculations.  Trajectory  prediction'! 
are  computed  in  the  main  program,  while  frequency-agility  calculations  and 
signal  processing  are  performed  in  subroutine  ERRO.  The  subroutine  GAIN 
calculates  antenna  pattern  data  (see  Section  E)  and  feeds  this  information 
to  ERRO.  The  remaining  subroutines  are  used  to  perform  statistical  analyses 
on  the  calculated  data  and  produce  output  plots. 

Various  versions  of  the  program  include  provision  for  different  beam¬ 
pointing  calculations  (including  linear  and  parabolic  predictions  and 
predictor-corrector  filters)  and  the  implementation  of  frequency  agility. 

The  trajectory  prediction  programs  fit  linear  and  parabolic  curves,  respec¬ 
tively,  to  a  set  of  three  data  points  to  predict  where  the  antenna  beam  should 
be  pointed  for  the  next  look.  The  frequency-agility  program  permits  selection 
of  start  and  stop  frequencies  and  the  frequency  step  size.  Each  frequency  is 
transmitted  for  three  pulses  and  then  stepped;  when  the  stop  frequency  is 
reached  a  new  cycle  is  begun.  Since  a  total  of  &8  pulses  are  transmitted 
for  each  look  at  the  target,  the  frequency  parameters  are  normally  chosen  such 
that  16  frequencies  are  transmitted  before  a  new  cycle  is  begun  in  order  to 
maximize  the  information  obtained.  Data  from  all  frequencies  may  be  averaged 
to  calculate  the  indicated  angular  position  of  the  target  or  adaptive-processing 
techniques  may  be  used  to  select  data  samples  from  the  available  set. 

F.  Results  of  Computer  Analysis 

From  the  computer  analysis  are  obtained  plots  of  the  elevation  tracking 
error  in  feet  as  a  function  of  time.  These  results  are  affected  by  antenna 
height  and  tilt  angle,  target  height  and  backscattering  characteristics,  radar 
frequency,  and  the  properties  of  the  reflecting  surface.  A  representative 
plot  of  tracking  error  is  shown  in  Figure  10;  this  particular  run  was  calculated 
for  an  isotropic-scatterer  target  having  a  radial  velocity  of  250  mph,  and 
was  initiated  at  a  distance  of  ten  miles.  The  frequency  of  operation  was 
chosen  to  be  5.5  GHz,  and  the  physical  face  of  the  antenna  was  vertical 
(antenna  normal  horizontal).  The  reflection  at  the  earth  surface  was  charac¬ 
terized  by  a  (voltage)  reflection  coefficient  of  0.5. 

The  beamshape  produced  by  the  antenna  varies  with  scan  angle,  beamwidth 
being  minimum  when  the  beam  is  normal  to  the  face  of  the  array,  and  broader 
for  non-zero  scan  angles.  Therefore,  the  physical  orientation  of  the  antenna 
may  be  chosen  so  as  to  minimize  antenna  beamwidth  and  maximize  tracking  accuracy 
for  the  region  of  most  interest.  Figure  11  shows  the  effect  of  tilting  the 
antenna  back,  which  enhances  accuracy  at  high  angles  at  the  expense  of  low- 
angle  accuracy  (compare  with  Figure  10).  The  most  desirable  tilt  angle  for 
the  antenna  depends  upon  a  number  of  factors  and  may  differ  with  mission 
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Figure  11.  Tracking  error  as  a  function  of  time.  Ha  =  100  ft,  11 
f  =  5.5  GHz,  p  =  0.5,  tilt  angle  =  30°. 


'W 


O  r 

requirements.  An  angle  of  30  was  chosen  arbitrarily  for  use  in  most  of  the 
remaining  computer  analyses. 

Figures  11  and  12  present  a  comparison  of  results  for  H  =  300  ft  and 
H  =  500  ft  for  an  antenna  tilt  angle  of  30°.  These  two  figures  illustrate 
the  influence  of  target  height  on  the  multipath-induced  tracking  errors.  The 
target  hfeight  is  seen  to  influence  both  the  amplitude  and  frequency  of  these 
multipath-induced  errors.  These  errors  are  also  affected  by  changes  in  antenna 
height:  Figure  13  shows  the  same  situation  as  Figure  11  (on  a  different  time 
scale),  but  with  antenna  height  reduced  to  15  feet.  The  peak  amplitude  of  the 
tracking  error  is  somewhat  less  at  the  lower  antenna  height,  and  the  frequency 
of  the  excursions  is  significantly  lower. 

Numerous  runs  of  the  type  shown  in  Figures  9-13  were  performed  in  order 
to  define  more  clearly  the  effect  of  multipath  returns  on  the  tracking  accuracy 
of  the  EAR  system.  An  antenna  tilt  angle  of  30  was  selected  for  most  of 
these  analyses,  and  an  isotropic-scatterer  target  having  a  radial  velocity  of 
250  mph  was  used.  A  surface  reflection  coefficient  of  0.5  was  chosen  for  all 
of  these  analyses.  Figure  14  summarizes  the  results  as  plots  of  peak-to- 
peak  tracking  error  as  a  function  of  target  height  for  two  values  of  antenna 
height  chosen  to  approximate  the  heights  of  the  initial  EAR  Test  Bed  site 
and  a  similar  antenna  mounted  on  a  tracked  vehicle.  Figure  14  indicates 
that  multipath  error  is  significant  for  targets  below  1200  feet,  and  is 
maximum  at  altitudes  of  about  400  to  800  feet. 

Figure  14  shows  that  tracking  error  due  to  multipath  can  be  of  such  a 
magnitude  as  to  severely  limit  radar  system  performance.  Various  methods 
for  alleviating  the  effect  of  multipath  returns  have  been  proposed,  and  are 
discussed  at  length  in  an  earlier  Georgia  Tech  report  [4],  The  most  practical 
method  for  implementation  with  the  EAR  is  frequency  agility,  and  the  appli¬ 
cability  of  this  approach  to  the  EAR  system  is  discussed  in  detail  in  the 
next  section. 
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Peak-to-peak  elevation  tracking  errors  as  functions  of  target 
height  for  15  and  100  feet  antenna  heights.  P  =  0.5  and 
frequency  is  5.5  GHz. 


IV.  APPLICABILITY  OF  FREQUENCY  AGILITY  TO  THE  EAR 


Frequency  agility  is  often  proposed  as  a  method  for  reducing  effects 
of  both  glint  and  multipath  on  radar  tracking  accuracy.  The  effectiveness 
of  frequency  agility  on  glint-induced  tracking  errors  has  been  investigated 
in  an  earlier  Georgia  Tech  report  [4].  While  the  EAR  presents  a  somewhat 
different  situation  (C-band  vs  X-Band),  the  results  set  forth  in  that 
report  are  generally  indicative  of  what  could  be  achieved  if  frequency 
agility  were  used  with  the  EAR  system  in  tracking  a  free-space  target. 

The  analysis  presented  in  Section  III  shows  that  the  effect  of 
multipath  returns  on  the  EAR  system  performance  is  significant  for  targets 
at  altitudes  less  than  1200  ft.  Numerous  methods  for  reducing  the  tracking 
error  have  been  proposed;  these  include  multiple-height  antennas,  space 
diversity  (in  the  azimuth  plane),  high-resolutioi.  antennas,  polarization 
agility,  shaped  antenna  beams,  and  frequency  agility.  The  discussions 
in  an  earlier  Georgia  Tech  Report  [5],  coupled  with  the  physical  constraints 
of  the  EAR  system,  make  frequency  agility  seem  to  be  a  promising  candidate 
technique  for  reducing  multipath-induced  tracking  errors.  Therefore,  an 
investigation  of  the  effectiveness  of  frequency  agility  in  reducing  these 
multipath-induced  tracking  errors  was  begun. 

When  frequency  agility  is  used,  a  set  of  pulses  is  radiated  in  sequence 
with  frequency  changed  from  pulse  to  pulse;  the  received  pulses  are  either 
averaged  to  obtain  a  more  accurate,  stable,  and  repeatable  indication  of 
target  position,  or  a  selection  rule  is  used  to  select  optimum  data  from  the 
available  set.  In  the  remainder  of  this  section,  the  term  "frequency  agility" 
will  be  used  to  describe  the  averaging  process,  and  the  term  "adaptive  pro¬ 
cessing"  used  to  denote  the  selection  process. 

Changing  the  transmitted  frequency  permits  the  acquisition  of  a  number 
of  independent  samples  of  target  location  by  changing  the  relative  phase 
between  direct  and  indirect  rays  (altering  their  electrical  path  lengths), 
by  changing  the  amplitude  and  phase  of  the  reflection  coefficient  of  the 
reflecting  surface,  or  by  affecting  the  scattering  characteristics  of  the 
target.  Within  the  normally  achievable  bandwidths  of  most  frequency-agile 
radars,  the  reflection  coefficient  of  many  reflecting  surfaces  remains 
essentially  unchanged.  The  effect  of  small  frequency  changes  on  target 
characteristics  has  not  been  exhaustively  investigated,  but  should  be  small 
for  most  low-altitude  targets,  particularly  in  the  vertical  plane  [5">.  Thus, 
the  principal  mechanism  by  which  frequency  changes  affect  tracking  error  is 
the  change  in  relative  phase  of  direct  and  indirect  signals  received  at  the 
antenna.  These  changes  in  electrical  path  length  produce  what  are  often 
substantial  shifts  in  the  location  in  range  of  the  extremes  of  the  tracking 
errors.  Figures  1?  through  19  show  the  tracking  errors  for  the  EAR  system 
with  operating  frequencies  of  5.3,  5.4,  5.5,  5.6,  and  5.7  GHz.  The  changes  in 
the  pattern  of  the  tracking  error  with  changes  in  frequency  are  clearly  seen, 
and  these  changes  indicate  that  practical  frequency  changes  may  indeed  be 
effective  in  obtaining  a  number  of  independent  samples  for  further  processing. 
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Figure  17.  Tracking  error  as  a  function  of  time. 

f  =  5.5  GHz,  p  =  0.5,  tilt  angle  =  30°. 
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Figure  18.  Tracking  error  as  a  function  of  time.  1^=  100  ft,  H 
f  =  5.6  GHz,  o  =  0.5,  tilt  angle  =  30  . 
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Figure  19.  Tracking  error  as  a  function  of  time.  H 
f  =  5.7  GHz,  p  ~  0.5,  tilt  angle  =  30°. 
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While  coherent  systems  such  as  the  EAR  may  usually  be  modified  for 
frequency  agility  with  a  minimum  of  transmitter  and  antenna  changes,  care 
must  be  taken  not  to  degrade  the  system  MTI  performance.  In  the  EAR  degrada¬ 
tion  is  circumvented  by  use  of  a  "triplet"  canceller  as  opposed  to  the  normal 
"sliding  window"  canceller  [6").  A  central  consideration  in  implementing  such 
a  system  is  the  amount  of  frequency  agility  bandwidth  required;  this  topic 
is  discussed  in  the  next  section. 


A. 


Required  Overall  Frequency  Agility  Bandwidth 
An  approximate  estimate  of  the  amount  of  overall  frequency  change 
required  may  be  arrived  at  by  considering  the  relative  phase  of  indirect  and 
direct  rays.  The  relative  phase.  9,  between  signals  which  traverse  two  paths 
whose  physical  lengths  differ  by  AR  is  given  by 


x  2tt  AR  _  „  f 

9  =  — - -  =  2tt  -  AR 

A  C 

When  frequency  agilty  is  used,  in  order  to  sample  all  possible  indicated 
target  positions,  it  is  necessary  to  change  9  by  n;  the  corresponding  required 
frequency  change,  Af,  may  be  arrived  at  as  follows: 


or 


9  +  TT  =  2rr— -f-  AR 
c 


Af  = 


2AR 


Taking  c  =  3  x  10  m/sec,  then 

^  v  in® 

Af  = 


2AR 

and  if  Af  is  MHz  and  AR  in  meters, 
Af  -  150/AR 


Since  for  multipath  situation  [7], 


AR  w 


2H  (H  +  II  ) 

a  a  t 
- 


the  required  frequency  change  becomes 

75f 


Af 


H  (H  +  H  ) 
a  a  t 


The  required  changes  for  II  =  15  feet  and  100  feet  have  been  calculated 
and  are  plotted  in  Figures  20  and  21.  These  data  indicate  that  the  400  to 
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Figure  20.  Frequency  change  required  Co  change  Che  relative  phase  of 
the  direct  and  reflected  signals  by  n.  Antenna  height  = 


Range  (meters) 


Figure  21.  Frequency  change  required  to  change  the  relative  phase  of 
the  direct  and  reflected  signals  by  rr.  Antenna  height  = 


500  MHz  potential  frequency  agility  bandwidth  of  the  EAR  system  is  sufficient 
to  be  at  least  partially  effective  in  reducing  multipath  effects  on  the  EAR, 
particularly  for  higher  altitude  targets,  higher  antenna  site,  and  shorter 
ranges . 

B.  Effect  of  Frequency  Agility  on  Angle  Tracking  Performance 

The  computer  analysis  of  EAR  tracking  performance  described  earlier  was 
used  as  tool  to  examine  the  effectiveness  of  frequency  agility  in  reducing 
multipath-induced  angle  tracking  errors.  The  format  selected  was  to  radiate 
16  frequencies  and  wait  0.1  seconds  between  looks  at  the  target.  The  indicated 
angular  errors  for  each  of  the  16  frequencies  were  then  averaged  together 
to  provide  an  estimate  of  true  target  position. 

The  degree  of  improvement  that  may  be  realized  with  frequency  agility 
is  a  function  of  target  and  antenna  heights,  and  of  the  frequency-agility 
bandwidth  (overall  frequency  excursion).  Data  from  a  number  of  runs  have 
been  summarized  by  presenting  peak-to-peak  tracking  error  as  a  function  of 
target  height  for  various  frequency-agility  bandwidths  and  antenna  heights. 

Figure  22  shows  one  such  set  of  data,  indicating  the  effectiveness  of 
various  frequency-agility  bandwidths  in  reducing  multipath-induced  tracking 
error.  The  three  curves  plotted  in  Figure  22  represent  the  performance 
achievable  with  fixed-frequency  operation,  200-MHz  bandwidth  frequency-agile 
operation,  and  400  MHz  bandwidth  frequency-agile  operation  for  a  100-foot 
antenna  height.  The  400-MHz  bandwidth  is  particularly  effective  in  reducing 
these  errors,  especially  for  higher  targets. 

As  discussed  earlier  and  shown  in  Figures  20  and  21,  the  antenna  height 
influences  the  amount  of  frequency  excursion  necessary.  Figure  23  shows  the 
results  obtained  when  frequency  agility  was  used  with  a  15-foot  antenna 
height.  While  tracking  performance  was  improved  with  frequency  agility,  the 
improvement  is  not  as  dramatic  as  seen  in  Figure  22  for  100-foot  antenna 
height. 

The  reflection  coefficient  also  influences  the  tracking  performance  of 
the  EAR  in  a  multipath  situation.  A  reflection  coefficient  of  P  =  0.5  was 
used  for  most  of  these  analyses,  however,  both  higher  and  lower  values  of 
p  were  also  used  during  the  course  of  the  analysis.  Figure  24  shows  the 
results  obtained  when  a  value  of  p  =  0.7  was  used,  indicating  the  while  the 
magnitude  of  the  errors  increases  with  increasing  p,  the  relative  effectiveness 
of  frequency  agility  remains  essentially  the  same. 

These  data  indicate  that  substantial  reduction  in  angle  tracking  error 
is  possible  when  frequency  agility  is  used  in  the  EAR  system,  provided  the 
full  frequency  agility  bandwidth  of  the  EAR  system  is  utilized. 

C.  Adaptive  Processing  for  the  EAR  System 

Use  of  the  relative  phase  difference  between  sum  and  difference  signals 
as  an  indication  of  the  acceptability  (or  quality)  of  the  tracking  information 
is  based  on  the  fact  that  when  tracking  a  single  target  in  free  space,  the 
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Figure  23.  Tracking  error  as  a  function  of  target  height  using 
frequency  agility.  H  =  15ft. 
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difference  signal  has  relative  phase  of  0  or  180° (depending  on  which  side  of 
the  axis  the  target  lies)  when  compared  with  the  sum  signal.  When  several 
unresolved  targets  (or  a  single  target  and  its  image)  are  illuminated  by  the 
radar,  this  phase  relationship  no  longer  holds.  The  basis  of  one  adaptive 
processing  scheme  is  to  use  this  relative  phase  between  sum  and  difference 
signals,  the  so-called  Complex  Indicated  Angle  (CIA)  [8l  as  a  measure  of  the 
quality  of  the  position  data.  A  plot  of  the  CIA  as  a  function  of  time  for  the 
same  run  shown  in  Figure  17  is  given  in  Figure  25.  A  correlation  is  seen 
between  the  phase  deviations  of  the  CIA  from  0  and  180  and  the  peaks  of  the 
tracking  error.  However,  these  phase  variations  are  rather  small  and  the 
presence  of  thermal  noise,  random  target  characteristics,  and  equipment 
inaccuracies  would  probably  make  it  difficult  to  appreciably  improve  the  quality 
of  the  track  data  by  using  the  CIA. 

Use  of  the  amplitude  of  the  received  sum  signal  as  a  means  of  identifying 
acceptable  track  data  has  also  been  proposed  [9].  Figure  26  shows  the  amplitude 
of  the  sum  signal  for  the  same  situation  shown  earlier  in  Figure  17.  Comparing 
these  two  figures  indicates  a  high  degree  of  correlation  between  the  minimum 
of  the  received  signals  and  the  peaks  of  large  tracking  error.  This  is 
probably  the  most  significant  correlation  which  may  be  used  to  assess  the 
quality  of  tracking  information,  and  on  which  has  been  proposed  to  reduce 
glint-induced  tracking  errors.  In  the  multipath  situation,  when  tracking  a 
complex  target,  variations  in  target  cross-section  with  frequency  may  be  as 
large  or  larger  than  the  variations  in  received  power  due  to  multipath  sienals. 
Thus,  there  is  no  guarantee  that  optimum  processing  is  obtained  by  selecting 
the  signal  of  minimum  amplitude.  However,  weighting  of  the  track  information 
based  on  the  amplitude  of  the  received  sum  signal  appears  to  be  a  reasonable 
approach,  since  it  has  a  high  likelihood  of  improving  the  overall  quality  of 
track  data,  if  for  no  other  reason  than  that  low-amplitude  signals  likely  to 
be  corrupted  by  thermal  noise  will  be  deemphasized.  Computer  analysis  of  more 
complex  target  models  in  a  multipath  situation  using  adaptive  processing  has 
shown  no  clearly  identifiable  improvement  in  the  quality  of  track  data  over 
conventional  frequency  agility.  While  there  were  cases  where  track  data 
improved  appreciably,  there  were  also  cases  where  it  did  not,  depending 
upon  details  of  the  lobing  of  the  target  and  the  lobing  due  to  multipath. 


(Degrees)  CIA  (Degree 


Signal  (dB) 


V.  IMPROVED  MTI  PROCESSORS 


A  number  of  factors  enter  into  che  design  and  specification  of  the  data 
processor  for  an  Mil  radar  system.  Among  the  constraints  are  requirements 
that  the  system  have  sufficient  clutter  attenuation  for  operation  in  a  heavy 
clutter  environment,  that  the  target  detectability  remain  relatively  constant 
for  the  range  of  expected  Doppler  frequencies,  and  that  processing  be  performed 
using  some  specified  number  of  received  pulses.  The  motivation  for  the  first 
two  requirements  are  rather  obvious,  and  the  last  requirement  is  dictated  by 
the  desire  to  minimize  the  number  of  pulses  required  on  a  given  target  (and 
consequently  to  maximize  the  number  of  targets  which  can  be  investigated) 
in  a  beam-agile  radar  such  as  a  phased-array  radar,  or  by  the  desire  to 
optimize  performance  when  frequency  agility  is  used  (by  minimizing  the  number 
of  pulses  transmitted  at  one  frequency). 

The  basic  points  to  be  developed  in  this  section  are  these:  (1)  The 
optimum  MTI  processor  (optimum  from  a  point  of  view  of  statistical  detection 
theory)  which  operates  on  more  than  two  received  pulses  has  not  yet  been 
developed.  (2)  If  such  a  processor  had  been  developed,  its  performance 
would  be  limited  by  practical  equipment  considerations.  (3)  Several  improved 
filter  designs  for  unstaggered  prf  systems  have  been  developed  which  offer 
substantial  clutter  attenuation  while  maintaining  more  nearly  uniform  response 
to  various  Doppler  frequencies.  (4)  Similar  improved  filter  designs  are 
needed  for  the  staggered  prf  case. 

In  this  section,  several  design  procedures  are  presented  for  the  realiza¬ 
tion  of  improved  MTI  processors  (cancellers,  or  clutter  filters),  and  charac¬ 
teristics  of  processors  designed  using  these  procedures  are  discussed.  Section 
A  reviews  receivers  that  are  optimum  from  the  point  of  view  of  statistical 
detection  theory.  Previous  work  directed  toward  development  of  optimum 
digital  MTI  processors  is  presented  in  Section  B,  where  difficulties  in 
obtaining  acceptable  MTI  processor  performance  are  discussed,  and  realistic 
constraints  set  by  equipment  limitations  and  by  system  performance  specifications 
are  outlined.  MTI  processors  that  are  optimum  from  a  filter-design  point  of 
view  are  developed  in  Section  G,  and  conventional  filter  design  procedures  are 
reviewed  in  Section  D.  In  Section  E,  several  design  procedures  for  use  in 
unstaggered  MTI  radar  systems  are  outlined.  Some  representative  results  of 
filter  responses  developed  using  these  design  procedures  will  be  presented, 
and  the  limitations  of  each  approach  discussed. 

MU  filters  for  use  in  systems  employing  staggered  prf  are  briefly 
discussed  in  Section  F.  The  effect  of  using  several  different  pulse  stagger 
ratios  ir.  a  staggered  prf  system  is  analyzed,  and  available  design  procedures 
for  staggered  systems  are  summarized. 


A.  Use  of  Statistical  Detection  Theory  in  Developing  Optimum  MTI  Receivers 
Development  of  optimum  receivers  has  been  of  substantial  interest  since 
radar  was  first  developed.  The  earliest  optimum  receiver  was  for  detection 
of  a  single  pulse  in  white  noi'-e.  This  concept  led  to  the  development  of  the 
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so-called  "matched  filter"  [10,11,12"),  namely,  one  which  maximizes  peak  signal- 
to-noise  ratio,  and  has  a  frequency  response  given  by  the  complex  conjugate 
of  the  voltage  spectrum  of  the  received  pulse. 

If  the  power  spectrum  of  the  noise  plus  received  clutter  varies  with 
frequency  (so-called  colored  noise),  then  the  optimum  filter  response  becomes 
(except  for  a  constant  time  delay)  the  complex  conjugate  of  the  voltage 
spectrum  of  the  received  pulse  divided  by  the  power  spectrum  of  the  receiver 
noise  plus  received  clutter  [13,14].  This  fact  was  used  by  Urkowitz  [ 15 ]  to 
derive  optimum  receivers  for  detection  of  targets  in  clutter. 

Rihaczek  [16]  has  pointed  out  that  the  class  of  filters  developed  by 
Urkowitz  is  optimum  only  when  thermal  noise  may  be  neglected,  and  that  the 
presence  of  both  fluctuating  clutter  and  thermal  noise  requries  more  complex 
filters  than  those  developed  by  Urkowitz. 

If  desired  targets  and  unwanted  clutter  returns  are  separated  in  time 
(range)  and/or  in  frequency,  substantial  improvement  in  performance  is  possible 
using  combined  signal  and  filter  optimization.  Delong  and  Hoffsteader  [ 17 , 18] 
consider  the  problem  of  detection  of  a  point  target  in  random  clutter  using 
combined  signal-receiver  optimization.  The  detection  of  a  target  of  known 
Doppler  shift  has  been  treated  by  Stuart  and  Westerfield  [19]  and  by  Van  Trees 
[20].  Spafford  [21,22],  Stutt  and  Spafford  [23],  and  Rummler  [24]  have  treated 
the  optimum  reciever  when  clutter  and  target  signals  have  different  areas  of 
occupancy  on  the  range-frequency  plane. 

In  many  cases  the  expected  Doppler  shift  of  the  received  signal  is  not 
known,  a  priori,  and  the  expected  range  of  signals  overlaps  the  clutter  in 
both  range  and  in  frequency.  The  optimum  estimation  receiver  for  this  case 
becomes  essentially  a  bank  of  matched  filters,  one  for  each  expected  Doppler 
frequency  [25];  this  configuration  is  very  similar  to  the  pulsed  Doppler  radar 
which  employs  a  comb  filter  or  a  filter  bank  followed  by  a  threshold  for  both 
velocity  estimation  and  target  detection. 

The  optimum  detection  receiver  corresponding  to  the  conventional  MTI 
radar  system  appears  to  have  been  first  discussed  in  the  radar  context  by 
Wainstein  and  Zubakov  [26].  Because  of  the  importance  of  this  work,  their 
basic  approach  to  this  problem  will  be  briefly  reviewed. 

The  formulation  of  the  optimum  WTI  receiver  is  one  of  testing  general 
Gaussian  hypotheses  for  the  case  of  a  nonfluctuating  target  and  interfering 
signals  which  are  Gaussian  random  variables.  Two  hypo  these  Hq  and  are 
defined  in  terms  of  the  received  signal  as  follows: 

Hq :  r (t)  =  n(t) 

r(t)  =  n(t)  +  m(t)  , 

where  m(t)  is  a  received  signal  reflected  from  a  point  target;  in  general 
m(t)  will  have  experienced  some  Doppler  shift.  The  interfering  signal  n(t)  is 
due  to  both  thermal  noise  and  reflections  from  clutter. 


Define  several  matrices: 


the  observation  matrix, 


the  mean  of  the  matrix  _r , 
m  =  E  ( r  )  , 

the  covariance  matrix 


A  -  Cov  (r)  =  E  £(r_  -  m)  ♦  (£T  -  mT )  J 


and  the  inverse  covariance  matrix 


a  "  a' 


The  optimum  Bayes  and  Neyman-Pearson  tests  are  both  likelihood 
ratios.  The  likelihood  ratio  (which  is  a  function  ■'f  both  the  Doppler  shift 
U)d  due  to  target  radial  velocity  and  the  initial  phase  of  the  received  signal  9) 
is  given  by  Van  Trees  [27]; 

lJ r (0,  wd)  =  exp  -  mT)  Q(r  -  m)  -  rTQrJ  J 

under  the  assumption  that  the  mean  value  of  n(t)  =  0. 

T  T  T 

Since  Q  is  symmetric  about  its  diagonal,  that  is,  Q  =  Q,  then  m  Qm  =  r  Qm. 

Using  this  fact,  the  likelihood  ratio  may  be  written 

L-jr  (0,  wd)  =  e^p  (mTQm)J  exp  (mT  Q  r) 

Since  both  the  anticipated  Doppler  shift,  u).,  and  the  initial  phase  8 
are  unknown,  the  average  likelihood  ratio  test  then  becomes,  assuming  all  u>d 
and  9  are  equally  probable, 

l|  (r)  =  f  (  exp  (mT  £  m)J  exp  (mT  £  r)  d^dG. 
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The  evaluation  of  this  likelihood  ratio  in  a  closed  form  is  a  formidable 
task.  The  first  integration  introduces  Bessel  functions  of  the  second  kind 
of  order  zero,  making  the  second  integration  difficult.  Wainstein  and  Zubakov 
[26]  have  applied  the  addition  formula  for  Bessel  functions  and  evaluated  this 
integral  exactly  for  the  case  where  the  observation  consists  of  two  pulses. 

The  optimum  MTI  receiver  derived  by  Wainstein  and  Zubakov  for  processing 
two  received  signals  consists  of  optimum  processing  of  both  the  in-phase  and 
quadrature  components  of  the  received  signal,  pairwise  subtraction  of  these 
two  in-phase  and  quadrature  samples,  formation  of  the  square  of  each  of  these 
differences,  and  comparison  of  the  sum  of  these  squares  with  the  threshold  [261. 
This  processing  corresponds  to  the  conventional  two-pulse  MTI  canceller. 

Selin  [28]  expanded  the  Bessel  function  of  the  second  kind  of  order  zero 
in  a  power  series  valid  for  small  ratios  of  signal  to  interfering  signal  in 
order  to  simplify  integration  of  the  Bessel  functions.  Selin  then  further 
confines  his  discussion  to  the  case  of  white  noise  interference  (which  is 
uncorrelated  from  pulse  to  pulse).  The  results  of  this  analysis  have  limited 
applicability  due  to  their  complexity  and  because  in  many  cases  of  interest, 
the  interfering  signal  is  highly  correlated  from  pulse  to  pulse  due  to  the 
presence  of  strong  clutter  returns. 

Brennon,  Reed,  and  Sollfrey  [29]  approximate  the  integral  of  the  Bessel 
function  by  a  finite  sum.  This  approximation  is  used  to  compare  the  perfor¬ 
mance  of  optimum  receivers  under  various  conditions,  but  it  does  not  give 
information  concerning  how  to  construct  an  optimum  MTI  receiver;  only  how  to 
approximate  its  performance  by  means  of  a  Doppler  filter  bank. 

From  this  review,  it  becomes  evident  that  the  specification  of  the 
optimum  MTI  processor  from  the  point  of  view  of  statistical  detection  theory 
is  a  formidable  task,  and  one  which  has  been  solved  exactly  only  for  the  case 
of  the  two-pulse  processor.  Because  of  the  difficulty  in  specifying  the 
performance  of  the  optimum  MTI  receiver,  considerable  attention  has  been 
focused  on  the  design  of  optimum  weighted  sums  for  processing  sampled 
sequences  of  the  return  from  moving  targets  in  a  clutter  environment.  This 
work  is  summarized  in  the  following  section. 

B.  Optimization  of  MTI  Processors  Using  Weighted  Sums  of  Sampled  Signals 

It  was  brought  out  in  the  previous  section  that  the  optimum  MTI  receiver, 
from  the  point  of  view  of  statistical  detection  theory,  is  known  only  for  the 
two-pulse  ease.  The  receiver  for • larger- numbers  of  received  pulses  has  often 
been  approximated  as  a  linear  combination  of  a  number  of  sample  values  of 
either  the  in-phase  or  the  quadrature  component  of  the  received  signal.  Maximi¬ 
zation  of  the  ratio  of  average  output  signal  (averaged  over  all  expected  values 
of  Doppler  frequency  shifts,  to.)  to  interfering  signal  has  been  treated  by 
Capon  [30].  Capon  shows  that  the  optimum  weight  functions,  defined  as  those 
which  optimize  the  average-output-signal-to-interfering-signal  ratio  (called 
the  reference  gain,  G  ,  or  more  commonly  the  MTI  improvement  I)  depend  only 
upon  the  covariance  matrix  of  the  interfering  signal.  For  highly  correlated 
pulse-to-pulse  interference,  such  as  that  due  to  slowiy  moving  clutter,  these 
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optimum  weight  functions  reduce  to  the  conventional  three-pulse  canceller  for 
the  case  of  processing  three  received  pulses.  Capon  also  shows  that  the 
reference  gain  for  the  three-pulse  canceller  closely  approximates  that  for 
a  large  number  of  received  pulses  when  processing  signals  in  a  background  of 
strongly  correlated  clutter. 

There  are  two  main  objections  to  Capon's  approach.  First,  there  is  no 
reason  to  believe  that  the  optimum  processor  may  be  realized  in  the  general 
configuration  assumed  by  Capon;  and  second,  the  concept  of  the  average  system 
gain,  G  ,  produces  poor  signal  detectability  at  some  Doppler  frequencies  of 
interest.  It  is  perhaps  appropriate  to  note  that  if  the  linear  processing 
format  discussed  by  Capon  were  the  configuration  of  the  optimum  processor, 
then  the  optimum  Neyman-Pearson  test  would  be  the  one  which  maximizes  G^ 

(see,  for  example,  Spafford  [22]). 

The  average  gain  G  ,  is  maximized  by  increasing  gain  at  frequencies  where 
clutter  return  is  small  and  decreasing  it  at  frequencies  where_clutter  is 
significantly  present.  Thus,  processors  designed  to  maximize  G  may  have 
unacceptably  low  responses  to  targets  with  Doppler  frequencies  m  the  same 
region  as  the  clutter.  The  problem  of  optimizing  the  response  of  MTI  systems 
for  a  wide  range  of  target  Doppler  frequencies  may  be  approached  by  considering 
the  processor  as  a  filter.  This  approach  is  developed  in  the  next  section. 

C.  MTI  Processor  Design  as  a  Filter  Optimization  Problem 

As  discussed  in  the  preceding  section,  maximization  of  the  average 
system  gain,  I,  is  often  not  a  very  satisfactory  method  for  optimizing  the 
processing  scheme  for  an  MTI  radar,  since  this  leads  to  poor  detectability 
of  targets  having  some  particular  range  of  Doppler  frequency  shifts.  This 
leads  one  to  consider  uniformity  of  response  of  the  processor  as  a  function 
of  tu  ,  as  an  important  consideration  in  system  design.  This  approach  leads 
naturally  to  considering  the  processor  to  be  a  filter  having  as  inputs  a  signal 
at  the  Doppler  frequency  and  a  signal  from  clutter  plus  thermal  noise  with 
known  power  spectral  density.  Then  the  filter  output  can  be  plotted  as  a 
function  of  input  Doppler  frequency;  one  such  representation  is  shown  in 
Figure  27  for  the  conventional  three-pulse  canceller.  As  one  can  see,  the 
response  is  very  non-uniform,  and  considerable  improvement  may  be  made  in  the 
shape  of  the  curve  while  still  maintaining  a  substantial  value  of  I. 

The  second  reason  why  maximizing  1  is  not  always  entirely  suitable  for 
optimizing  system  response  is  that  there  are  often  substantial  practical 
limitations  which  are  not  included  in  the  theory.  Most  modern  high-performance 
MTI  processors  utilize  digital  processing,  mainly  to  obtain  storage  without 
recirculating  analog  delay  lines.  While  round-off  error  in  the  quantization 
or  digitizing  process  does  not  usually  appreciably  limit  performance,  analog- 
to-digital  quantization  errors  usually  constrain  I  to  be  substantially  less 
than  its  theoretically  achievable  maximum  value.  A  comprehensive  treatment  of 
the  limitations  in  improvement  due  to  quantization  errors  is  given  in  Appendix  F. 
Therefore,  it  is  often  possible  for  I  to  be  reduced  somewhat,  providing  more 
uniform  detectability  of  targets  of  different  Doppler  shifts,  without  substan¬ 
tially  affecting  overall  system  performance. 
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filter  as  a  function  of  Doppler  frequency. 


Another  factor  which  limits  the  performance  of  some  real  radar  systems 
is  the  presence  of  slowly  moving  discrete  targets  such  as  birds,  insects,  and 
automobiles  [31,32].  The  presence  of  these  extraneous  targets  sometimes 
requires  that  a  "stop-band"  be  established,  usually  centered  about  aero 
frequency,  in  order  to  reject  these  unwanted  returns  which  would  otherwise 
completely  overload  the  system. 

Therefore,  one  might  conclude  from  these  remarks  that  a  practical  design 
procedure  for  digital  MTI  filters  would  be:  (1)  establish  a  desired  and 
achievable  value  of  I,  based  both  on  desired  system  performance  and  on  practical 
equipment  limitations,  (2)  within  this  constraint,  produce  an  optimally  uniform 
response  for  all  Doppler  frequencies  of  interest,  while  (3)  rejecting  unwanted 
returns.  The  filter  responses  are  a  function  of  the  shape  of  the  clutter 
spectrum,  the  optimization  criteria  applied  to  the  response  (such  as  minimum 
rms  error,  equal  ripple  in  the  pass  band,  etc.)  and  the  expected  range  of 
velocities  of  the  desired  and  undesired  targets. 

The  following  discussions  will  be  largely  confined  to  consideration  of 
transversal  filters  [33,34]  (nonrecursive  filters  or  those  with  no  internal 
feedback  loops),  because  in  a  real  application  only  limited  number  of  pulses 
may  be  processed  from  each  target.  Several  constraints  determine  the  number 
of  pulses  that  may  be  processed  from  a  given  target.  In  a  beam-agile  radar 
such  as  a  phascd-array  system,  minimizing  the  number  of  pulses  on  a  given 
target  maximizes  the  number  of  targets  the  radar  can  accomodate.  If  frequency 
agility  is  used,  the  radar  must  remain  at  a  given  frequency  for  a  sufficient 
number  of  pulses  to  extract  the  desired  information  concerning  a  target; 
minimizing  the  number  of  pulses  on  target  thus  maximizes  the  number  of 
available  frequencies  the  radar  may  radiate  in  a  specified  time.  The  perfor¬ 
mance  of  several  conventional  transversal  digital  MTI  filters  is  discussed 
as  Appendix  A.  While  a  recursive  filter  (one  containing  feedback  loops)  could 
be  used  and  its  transient  response  truncated  after  the  desired  number  of  pulses, 
the  response  of  such  a  truncated  recursive  filter  may  always  be  realized  as 
a  transversal  filter.  The  difference  between  the  two  lies  in  the  practical 
implementation  of  the  filter. 

D.  Conventional  Digital-Filter  Design  Procedures 

Various  methods  have  been  developed  for  designing  digital  filters  [35,3b"1. 
Their  design  is  often  approached  by  defining  an  analog  filter  prototype  and 
appropriately  transforming  the  response  to  obtain  the  z-transform  of  the  desired 
filter.  In  general,  this  approach  yields  a  recursive  filter;  while  this 
filter's  output  may  be  truncated  after  the  desired  number  of  pulses,  there  is 
generally  little  control  over  the  number  of  pulses  required  to  closely  appro¬ 
ximate  the  desired  steady-state  response. 

To  illustrate  the  errors  in  filter  response  that  may  occur  due  to 
truncation  of  a  recursive  filter  designed  for  a  certain  steady-state  response, 
a  four-pole  Butterworth_ filter  response  was  considered.  Its  z-transform  was 
expanded  to  powers  of  z  1  by  long  division  and  the  series  truncated  after 
a  selected  number  oi  terms.  The  impulse  response  of  the  filter  represented  by 
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this  series  was  then  calculated  to  determine  the  truncated  frequency 
response.  Figure  28  compares  the  steady-state  response  and  responses 
obtained  by  truncating  the  filter  response  after  3,  5,  and  7  pulses.  As 
can  be  seen,  the  response  of  the  truncated  series  is  a  poor  extreme  case, 
due  to  the  rapid  low-frequency  roll-off  of  the  filter,  but  it  serves  to 
illustrate  the  need  for  specialized  design  procedures  where  the  number  of 
available  samples  is  limited. 

Three  procedures  are  commonly  used  in  designing  transversal  digital 
filters  [35-43].  The  first  involves  specifying  the  sampled  finite 
impulse  response  of  a  filter  (obtained  as  an  input  specification  or  as  a 
transform  of  a  frequency-sampled  response)  and  utilizing  an  impulse-invariant 
transformation  to  specify  the  digital  filter.  Various  weighing  windows  may 
be  used  to  smooth  the  ripples  in  the  resulting  frequency  response.  The  second 
method  uses  a  Fourier-series  approximation  (with  appropriate  windows  often 
incorporated)  to  a  desired  response,  while  the  third  method  is  a  direct-search 
method  using  linear  programming  techniques  to  optimize  the  desired  response. 
These  techniques  are  not  particularly  applicable  to  the  design  of  MTI  processing 
filters,  since  they  do  not  incorporate  clutter  characteristics  into  the  filter 
design  procedure  as  a  design  specification.  In  many  cases  they  accept  a 
specific  filter  shape  as  the  design  goal,  rather  than  developing  optimum 
processing  for  a  predetermined  number  of  received  pulses,  a  constraint  which 
is  determined  by  the  energy  budget  of  the  overall  radar  system.  Four  techniques 
that  have  been  developed  for  design  of  optimum  MTI  processing  filters  for 
radars  using  unstaggered  prf's  will  be  outlined  in  the  next  section.  The 
staggered  prf  case  is  considerably  more  complex;  however,  some  related  work 
has  been  done  in  this  area  and  will  be  reviewed  in  Section  F.- 

E.  Optimum  Transversal  Digital-Filter  Design  Procedures  for  Unstaggered  prf 
Systems 

The  design  of  optimum  MTI  filters  should  take  into  account  a  number  of 
factors,  including  maintaining  desired  values  of  I,  providing  uniformity  of 
response  for  some  range  of  Doppler  frequencies,  and  in  some  cases  providing 
a  stop-band  (usually  centered  about  zero  frequency)  for  rejection  of  unwanted 
targets.  In  this  section  four  different  design  procedures  are  described  which 
take  some  of  these  factors  into  account. 

1.  "Cost"  Minimization  MTI  Filters 

Jacomini  (.44]  has  developed  a  design  procedure  that  takes  into  account 
the  clutter  attenuation  in  the  pass  band,  the  response  in  a  stop-band,  and 
ripple  in  the  pass  band.  A  "cost"  is  assigned  to  each  of  these  factors  and 
designs  developed  which  minimize  the  overall  "cost"  of  the  filter.  This 
approach  has  the  disadvantage  that  there  is  usually  no  straightforward  means 
for  selecting  the  various  costs.  In  practice,  when  the  filter  is  being 
designed,  the  improvement,  I,  is  constrained  to  be  some  constant,  the  cost 
is  apportioned  between  the  pass  and  stop  bands,  and  the  relative  costs 
adjusted  until  an  intuitively  satisfying  filter  response  is  obtained. 
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More  straightforward  means  for  obtaining  the  same  types  of  results  are  described 
in  following  sections. 

2.  Filter  Design  by  Linear  Programming 

The  frequency  response  of  transversal  filters  can  be  completely  specified 
by  an  equation  which  is  a  linear  function  of  the  C's  (see  Appendix  F) .  This 
fact  simplifies  the  design  of  this  class  of  filters,  permitting  use  of  well- 
established  linear  programming  techniques.  Optimization  by  linear  programming 
makes  use  of  the  fact  that,  given  a  set  of  linear  constraint  equations,  along 
with  an  "objective  function"  which  is  also  a  linear  equation,  the  extreme 
value  of  the  linear  objective  function  lies  on  the  boundary  of  the  "feasible 
region"  which  is  defined  by  this  set  of  constraint  equations.  This  important 
fact  gives  rise  to  the  "simplex"  method  of  linear  programming;  an  optimization 
scheme  which  methodically  searches  the  boundary  of  the  feasible  region  for 
the  extremum  of  the  linear  objective  function. 

When  a  desired  filter  response  is  to  be  approximated,  several  techniques 
may  be  used  to  determine  the  appropriate  filter  coefficients:  two  of  these 
are  the  "least-mean-square-error"  (LMSE)  method  and  the  "minimax"  method. 

The  LMSE  minimizes  the  least-mean-square  error  between  the  specified  response 
and  the  approximation  for  all  frequencies  of  interest. 

In  the  LSME  method,  if  Q(m)  is  the  specified  response  function  and  G(u)) 
is  the  response  of  the  filter,  the  goal  is  to  minimize 

\  [G(<d)  -  Q(<o)]2  dou  , 
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where  R  is  the  frequency  range  of  interest  [451.  For  many  applications  this 
method  of  filter  design  is  desirable  and  is  treated  in  Section  3.  However, 
simplex  methods  are  not  applicable  for  LSME  specification  since  the  resulting 
equations  are  not  linear. 

The  "minimax"  method  is  used  to  reduce  the  maximum  deviation  from  the 
filter  specification  at  a  number  of  selected  frequencies.  The  minimax  formula 
reduces 


for  a  number  of  discrete  values  of  u)  contained  in  R,  the  frequency  range  of 
interest.  For  example,  a  typical  application  might  be  to  design  MTI  filters 
with  bounded  ripple  over  some  specified  frequency  range  while  maximizing 
attenuation  for  undesired  targets  and  clutter. 

Such  a  minimax  program  was  implement.ed  on  the  Univac  1108  using  a  well- 
known  linear  programming  scheme  [46]  The  procedure  used  was  to  identify 
a  number  of  frequencies  at  which  the  frequency  response  would  be  controlled. 

At  each  of  these  frequencies,  minimum  and  maximum  values  of  frequency  response 
were  specified.  In  addition,  a  linear  objective  function  (LOF)  was  defined. 


For  this  work,  the  clutter  attentuation  divided  by  C  was  chosen  as  the  LOF; 
since  ~  1  for  most  practical  filters,  minimizing  ?lie  LOF  corresponded  to 
maximizing  the  clutter  cancellation.  The  linear  programming  technique  which  was 
used  minimized  the  LOF,  subject  to  the  constraint  that  the  frequency  response 
remain  within  the  bounds  established  earlier.  One  substantial  problem 
encountered  using  this  technique  was  that  roundoff  error  seriously  influenced 
the  clutter  attenuation  that  could  be  calculated.  Nevertheless,  practical 
filter  responses  could  be  developed  using  these  techniques.  A  second 
difficulty  involved  the  large  number  of  equations  that  had  to  be  entered  into 
the  program.  Successful  results  were  obtained  for  three-  and  five-pulse 
filters,  but  greater  than  five  pulses  was  beyond  the  scope  of  this  program. 

This  linear  programming  scheme  has  been  used  to  design  a  number  of  filter 
responses.  Most  of  the  results  discussed  here  will  be  for  five-pulse  filters. 
The  filter  response  for  these  five-pulse  filters  were  specified  in  250-Hz 
intervals  from  zero  to  2500  Hz,  and  the  constraint  at  zero  frequency  was  that 
the  response  be  non-negative. 

Figure  29  shows  three  of  the  MTI  filter  responses  were  developed  using 
the  program.  The  tradeoff  between  ripple  and  low-frequency  clutter  attentuation 
can  be  seen  clearly.  The  3-dB  filter  has  sharper  low-frequency  cutoff  than 
the  0.8-  and  1.8-dB  filters,  thus  increasing  the  attenuation  for  clutter  and 
slowly  moving  targets.  With  more  available  pulses,  the  cutoff  slope  could 
be  further  increased  (thus  increasing  clutter  attenuation)  while  retaining 
constant  ripple. 

The  ripple  constraint  in  the  pass  band  can  be  made  so  restrictive  that 
the  filter  provides  essentially  no  clutter  attenuation,  as  is  seen  in  Figure 
30.  When  the  ripple  constraint  for  this  particular  filter  was  changed  from 
3  dB  to  0.8  dB,  it  failed  to  provide  appreciable  attenuation  for  low  frequencies 
The  minimum  specified  frequency  was  250  Hz. 

This  relationship  between  ripple  and  clutter  attenuation  was  investigated 
for  five-pulse  radar  filters,  resulting  in  the  graph  shoxfn  in  Figure  31.  The 
curves  show  the  ripple  vs  clutter  attenuation  relationship  for  500  Hz  and 
750  Hz  filters  calculated  using  a  36-bit  computer.  Extensions  to  these  curves 
and  slightly  increased  accuracy  could  have  resulted  if  more  computing  accuracy 
had  been  available.  It  is  possible  to  project  the  clutter  attentuation  for  a 
wide  variety  of  filters  by  using  curves  of  this  type. 

All  the  filters  previously  discussed  have  been  developed  to  suppress 
clutter  having  a  spectral  width  that  is  narrow  in  comparison  with  the  system 
bandpass;  the  clutter  has  been  treated  essentially  as  ar  impulse  at  zero 
frequency.  However,  situations  often  arise  in  which  targets  or  clutter  lying 
somewhat  higher  in  frequency  need  to  be  suppressed.  This  may  be  caused  by 
unusual  clutter  conditions  or  by  the  presence  of  undesired  slowly  moving  targets 
A  stop  band  which  will  suppress  these  frequencies  can  be  synthesized  easily 
by  simply  rewriting  the  system  constraint  equations  to  define  these  new 
requirements.  This  will  usually  require  that  the  ripple  constraint  be  relaxed 
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Cutoff  =  500  Hz 


Figure  31.  Ripple  vs.  clutter  attenuation  for  5-pulse  MTI  filters  with  500  Hz 
and  750Hz  lower  cutoff  frequencies. 


somewhat.  In  Figure  32  some  low-frequency  stopband  filters  with  ripples  of 
1,  2,  and  3  dB  which  were  developed  using  linear  programming  techniques  are 
plotted.  Not  shown,  however,  are  the  low-frequency  attenuation  characteristics 
below  30-dB  attenuation.  The  3-dB  ripple  filter  achieved  approximately  10  dB 
more  attenuation  in  the  stop  band  than  the  2-dB  ripple  filter.  Similar 
differences  were  noted  between  the  1-  and  2-dB  cases. 

These  cases  show  the  relationship  which  exists  between  ripple  and  clutter 
attentuation  for  some  practical  filters,  give  guidelines  for  filter  performance, 
and  demonstrate  the  use  of  filters  with  stop  bands.  In  all  of  these  filters, 
performance  was  ultimately  limited  by  the  number  of  pulses  available,  end 
processing  larger  numbers  of  pulses  would  have  resulted  in  filters  having 
more  desirable  performance. 

3.  MTI  Clutter  Rejection  Filters  Using  an  rms  Error  Specification 

The  specific  problem  investigated  was  to  maximize  detectability  of 
moving  targets  over  some  range  of  Doppler  frequencies  while  maintaining  some 
minimum  clutter  attenuation.  The  procedure  used  was  to  constrain  the  target- 
to-clutter  ratio  improvement  factor,  I,  to  be  some  selected  value  while 
minimizing  the  function 

rV2  _ 2 

(  r  (G(f>  -  F (f);  df  , 

where 

f  =  the  pulse  reprtition  frequency, 

D  =  a  fraction  0<f|<0.5, 

G(f)  =  the  power  response  of  the  filter,  and 


G(f)  =  average  response  of  the  filter. 

The  choice  of  T)  amounts  to  establishing  a  pass  band  of  interest;  for  the  EAR 
parameters,  1]  =  0.1  corresponds  to  minimizing  the  error  for  targets  above 
approximately  30  mph.  The  value  of  G(f)  was  chosen  equal  to  one. 

The  problem  now  posed  is  one  which  may  be  solved  readily  by  Lagrange 
multiplier  techniques.  The  steps  involved  in  solution  of  the  problem  are 
(1)  write  an  expression  for  the  improvement,  I,  (2)  write  an  expression  for 
the  error,  (3)  apply  Lagrange  multiplier  techniques,  and  (4)  solve  for  the 
coefficients.  These  steps  are  outlined  in  the  following  sections  It  is 
worth  noting  a  related  analysis  has  been  performed  by  Martin  [33],  for  cases 
considered  here,  his  analysis  is  a  special  case  of  the  work  which  follows  with 
T]  =  0. 
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250  Hz 


a.  Improvement  for  an  Arbitrary  N-Pulse  Canceller 
As  shown  in  Appendix  F,  the  power  response  is 

N 

GGa)  =  c  +2  '2L  c  cos 

o  q 

q=i 

a..d  a  Gaussian  clutter  power  spectral  density  input  is  assumed, 


2  2 

„  _  c  -r/20 

CT  = -  e  c 

c  /2tt 
c 

The  clutter  output  is  given  by  CouC,  where 


f  / 2  N  /  p  -f^/2rr^\ 

Cout  =  L  (Co  +  2  ^  Cq  C0S  qu,T)  W7Tr.  e  c)df 


out  '  0  *rr.  q 

-f  /2  q-1  n  c 

r 


£  /2  N  ,  r  P2/?  2\ 

=  2;r  (C  +  2  Z. ,C  C0S  q»T>  (jyg;  e  c)  df  . 
o  q=l  c 


Since  ac«fr> 


cw  LN 

c  «  f  (c  +  2  JT  c 

out  J  0 

o  q-1 


,  e'£2/2o")«  ■ 


cos 


or 


=out  -  jk  K  v-£2/2°= dE + 2  4  f  %  ~ 

C  O  vj— ■!  ^ 


r-  N  „  2  2  2  2-, 

C  =  C  fc  +  2  Z  c  e  q  "c  ■ 
out  o  q  _ 


which  is  the  desired  result. 

The  average  target  output  power,  T  ,  is 


T.  f 

T  =  f  r  G(f)  d£ 

O  r 

r  o 

T.  t_  p  N 


=  ~ir  C  +  2  r  C  cos  2rrfqT 
£  in  ,  q  _ 

q=l  ^ 


£  J  i_  o 
r  o 


df 


-  C  T. 
0  l 
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Then 


I  = 


N  0  2  2  2  2 

C  +  2  f  C  e“2q  n  T  °c 
0  q=i  q 


This  expression  checks  with  earlier  results  if  the  exponentials  are 
expanded  in  a  Taylor  series  and  the  first  non-zero  terms  retained.  Consider 
the  three-pulse  canceller  where  Cq  =  6,  =  -4,  0,^  =  1.  Then 

.  ir  0  2m2  2  q  22  2n 

1  _  If,  q  -2n  T  a  ,  0  -8tt  T  a  j 

T  -  6| .  6  •  8  e  c  +  2e  CJ 


X  X 

substituting  e  «  1  +  X  +  -r-  , 


I  = 


o  4  4 
Sn  ac 


as  before. 


If 


f  =  5  x  10 
r 

O  =0.22 
r 

X  =  5.5  cm 

are  used  to  calculate  I  for  the  three-pulse  canceller  (Cq  =  6,  =  -4, 

and  =  1) >  a  value  I  =  83.01  dB  is  obtained,  which  agrees  closely  with 
earlier  results. 

b.  Error  Statement 

The  error  we  have  chosen  to  minimize  over  the  reg?cn  Tlf  to  f  /2 

r  r 

is  the  function 

f  /2  , 

E  =  J  r  (G(f )  -  1)  df  , 

Tlf 

r 

Now, 


[G (o>)  -  l]2  =  [C  -1  +  2  r  C  cos  qojT]' 


q=l 


60 


=  (C  -  1)  +  4C 

0  o  r  q 

q=l 


N  N 

C  cos  qajT  -  4  £2  cos  qu>t  + 

q=l 


N  N 

*z  r  C  C  cos  qojT  cos  kwT 
q=l  k=l  q  R 


Therefore, 


f 


E  =  f  r  [G(f)  -  IT  df 
Ttf 


-  (C0  -  l)2(fr/2)(l  -  2D)  -2C0  X  J:  c  y 


q=l  k=l 


q  k  r  qlc 


N  N  N 

2  H  _fl  r_  sin  2i™<1  +  4  ZL  CnCkfraak 

q-1  nq  q-1  k-1  q  R  *  qK 


where 


O' 


qk 


C  sin  2TT'P(g  -  k)  _  sin  2nT)(g  +  k)  j  . 

I  "  4rr(q-k)  4n(q  +  k)  q  r 


if,  07>  sin  4rfc'f,_1 
_^1  ..  2T)  -  —  -  .  ^ 


q  =  k 


c.  Optimization 

Using  stands  d  Lagrange  multiplier  techniques,  we  wish  to  minimize 


F  =  E  ■  'f  ■  v)  f  C 
VI  /  r  o 


subject  to  the  constraint  that 
j-  =  constant  =  Y 

This  is  done  by  solving  the  system  of  equations 
dF 


9C 


=  0 


5F 

ac. 


=  o 


(f-v). 


bl 


Noting  that  or  .  -  O'  ,  these  expressions  are  somewhat  simplified;  for  a 
fhree-pulse  filter,  this  becomes 

Cq(1  -  27])  -  +  x(l  -  v)  =  1  -  27] 

«iin  2-rTI  o_ZmZ  Z.  —  2  sin  2ttT| 

-2C  +  c  8a  +  C08a10  +  X(2e"2n  T  V  "2  rr 

OTT  111  /.  1 2. 


2  2  2 


■“.TT1  +  Cl8°21  +  C2to22  +  X(2e'8"  T  °c>  “  ‘2  ^^T21 


„  22  2  0  2_2  2 
(1  -  y)  Cq  +  C1  2e"2rI  T  °c  +  C2  2e"8TT  T  ac  =  0 


The  system  of  equations  readily  expands  for  larger  numbers  of  received 
pulses.  A  computer  program  was  written  to  solve  for  the  C's  for  up  to  100 
pulses.  Results  obtained  using  this  program  are  discussed  in  the  next  section. 

d.  Results 

The  filter  responses  for  the  conventional  three-pulse  MTI  filter 
and  one  of  the  minimum-error  filters  designed  for  I  =  60  dB  is  shown  in 
Figure  33.  The  improvement  in  detectability  for  targets  having  low-frequency 
Doppler  returns  is  evident. 

The  shape  of  the  filter  response  is  a  function  of  the  exact  value  of  I 
specified;  Figure  34  shows  results  for  optimum  . asponses  for  I  =  10  dB, 

30  dB,  and  60  dB,  indicating  the  increased  ripple  associated  with  larger 
values  of  I.  This  illustrates  the  desirability  of  choosing  I  no  greater 
than  necessary  because  of  the  consequent  compromises  in  filter  responses. 

The  flatness  of  the  filter  increases  as  the  number  of  pulses  processed 
increases,  but  the  number  of  ripples  also  increases.  This  is  illustrated  in 
Figure  35  which  shows  responses  for  I  =  60  dB  for  3,  5,  and  7  pulses  processed. 

Differences  in  responses  for  Tj  =  0  and  7]  =  0.1  are  not  particularly 
great,  but  the  case  7]  =  0.1  does  show  a  slight  improvement  in  flatness  of 
response.  Figure  36  compares  results  for  5-pulse  response  with  1  =  60  dB 
for  T]  =  0  and  7]  =  0.1. 
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Figure  33.  Comparison  of  the  response  of  a  conventional  three-pulse  MTI 
filter  and  a  minimum  rms  error  filter  designed  for  I  =  60  dB, 


34.  Comparison  of  responses  of  minimum  rms  filters  designed  for 
I  =  10,  30,  and  60  dB.  T)  =  0.1,  prf  =  5000  pps ,  and  crc  =  8 
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filters  for  7}  =  0  and  T>  =  0.1.  I  =  60  dB,  a  =  8.0,  prf  =  5000  pps 


4.  Maximally  Flat  Non-Recursive  Digital  MTI  Filters 

The  digital  MTI  filters  designed  up  to  this  point  have  had  ripples 
in  the  pass  band,  with  consequent  variation  in  target  detectability.  This 
variation  is  objectionable  in  some  cases;  therefore,  maximally  flat  non¬ 
recursive  digital  filter  responses  having  a  specified  improvement  factor  have 
been  developed  for  consideration  for  implementation  in  the  EAR  radar  system. 

By  a  "maximally  flat"  filter  is  meant  one  having  a  number  of  its  deriva¬ 
tives  with  respect  to  frequency  set  equal  to,„zero  at  some  specified  frequency 
or  frequencies  [4l].  For  the  MTI  case,  it  is  also  required  that  these  filters 
have  a  certain  specified  amount  of  improvement  in  target-to-clutter  ratio. 

In  addition,  the  filter  must  have  a  non-zero  amplitude  specified  at  a  given 
frequency,  in  order  to  prevent  a  solution  which  is  identically  equal  to  zero. 

The  procedure  which  was  used  to  derive  the  maximally  flat  nonrecursive 
digital  MTI  filters  is  described  below.  As  before,  the  power  response  of  the 
filter  is  given  by 

N 

G(cu)  =  C  +  2  £  C  cos  qu>T  , 
o  q=l  q 

and  the  target-to-clutter  improvement,  I,  is  given  by 


+  2l 


_  ,  „  2  22  2. 
C  exp  (-2q  n  T  a  ) 
q  c 


Following  the  above  discussion,  choose  a  value  for  I  =  — ,  and  choose  GCnf^)  =  1 
Then  there  are  N-2  derivatives  which  may  be  set  equal  to  zero.  The  points 
at  which  to  set  these  derivatives  equal  to  zero  were  chosen  to  be  to  =  0  and 
(U  =  rtf  .  The  way  that  these  derivatives  are  apportioned  between  these  two 
frequencies  determines  the  shape  of  the  filter  characteristic. 

A  three-pulse  canceller  has  only  one  derivative  to  set  equal  to  zero;  in 
order  to  obtain  a  high-pass  filter  characteristic  this  derivative  must  be  zero 
for  (jd  =  nf  *  For  the  three-pulse  canceller,  the  equations  for  the  C's  become 

G(TTfr)  =  1 


=  U)  “  TTf 


The  filter  response  for  I  -  60  dB,  a  =8,  and  T  =  1/f  =  1/5000  is  shown 

in  Figure  37.  Figure  38  compares  this  maximally  flat  three-pulse  filter 
response  with  that  of  a  minimum-rms-error  three-pulse  canceller,  indicating 
the  improved  flatness  of  response  obtained  using  a  maximally  flat  filter  response. 
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If  N  pulses  are  processed,  there  are  N-2  derivatives  which  may  be  set 
equal  to  zero,  and  these  derivatives  are  specified  at  m  =  0  and  tjt>  =  TTf  . 

Define  the  order,  i,  of  the  filter  as  being  the  number  of  derivatives  set 
equal  to  zero  at  u>  =  0.  Responses  for  filters  with  i  =  0,  and  N  =  3  and  N  =  4 
are  compared  in  Figure  39,  where  it  may  be  seen  that  increasing  N  has  little 
effect  on  the  filter  response.  However,  higher-order  filters  for  N  >  3, 
permit  establishment  of  a  stop  band  to  reject  slowly  moving  targets.  Responses 
for  N  =  4  and  N  =  5  for  various,  values  of  i  are  shown  in  Figures  40  and  41  to 
illustrate  the  shape  of  the  pass  and  stop  bands  produced  by  various  choices 
of  N  and  i. 

F.  MTI  Processors  Using  Staggered  prf 

One  of  the  serious  limitations  of  the  systems  discussed  up  to  this  point 
is  that  the  filter  responses  exhibit  so-called  "blind  speeds"  --  target 
radial  velocities  for  which  the  filter  output  is  zero.  These  speeds  are  those 
for  which  the  Doppler  frequency  is  an  integral  multiple  of  the  system  prf. 

If  it  is  appropriately  implemented,  frequency  agility  may  be  used  to 
reduce  effects  of  blind  speeds  on  system  performance.  An  analysis  of  the 
effects  of  frequency  agility  on  blind  speeds  is  presented  as  Appendix  D,  where 
it  is  shown  that  this  technique  is  effective  only  against  targets  having  large 
radial  velocities. 

The  other  method  for  reducing  the  effects  of  blind  speeds  is  the  use  of 
prf  stagger,  that  is,  varying  the  interpulse  spacing  from  interval  to  interval. 
Use  of  stagger  effects  both  the  shape  of  the  frequency  response  and  the  MTI 
improvement,  I,  of  the  system.  An  analysis  of  effects  of  prf  stagger  on 
conventional  MTI  processors  will  now  be  briefly  discussed. 

Changing  the  interpulse  period  by  a  fraction,  e,  will  greatly  reduce  the 
depth  of  the  nulls  of  the  filter  response  and  provide,  in  some  cases,  a 
desirable  response  without  blind  speeds. 

For  a  three-pulse  filter  having  a  stagger,  e,  about  the  center  pulse 
the  impulse  response,  h(t),  is  given  by 

h(t)  =  6 (t  +  T(l-e))  -  2  6(t)  +  6(t  =  T(l+e)) 

and  the  frequency  response,  G(uj),  of  this  filter  is  given  by  the  Fourier 
Transform  of  h(t), 

G(m)  =6-8  cos  u)T  cos  ojTe  +  2  cos  2u)T 

Figure  42  shows  the  response  of  such  a  filter  with  values  of  e  ranging 
from  0  to  0.1,  showing  how  the  lower  blind  speeds  are  eliminated. 

The  reduction  in  I  due  to  staggered  prf  may  be  determined  as  shown  by 
Nathanson  [48]. 


3  pulse 


Figure  41.  Maximally  flat  filter  responses  for  a  five-pulse  MTI  filter 
for  i  =■  0,  1.  and  2. 


The  residue,  E^vt),  from  a  conventional  three-pulse  canceller  is 
Er(t)  =  E(t)  -  2E(t-T)  +  E(t-T1-T2) 

where 

E(t)  =  sample  at  time,  t 
T  =  T (1  +  e) 

T2  =  T ( 1  -  -:) 

Then  the  mean  square  value  is  given  by 

E^(t)  =  [E(t)  -  ZECt-T^  +  E(t-T1-T2]  , 

which  reduces  to 

=  6E2(t)  -  |-CR(T1)  +  R(T2)]  +  i  R(T1  +  T2)  . 

Let  the  normalized  autocorrelation  be  defined  as 


Therefore , 


The  improvement,  Ig,  with  staggered  prf  is  I  multiplied  by  the  average 
gain  (six) 


S  1  -  3CP (Tx)  +  PCKT^)  +  jpCTjCl  +  K)] 

The  loss  due  to  staggering  is  [40] 

I  -  4aT  a  -  K)2  16  [k2  - 

-r  ~  f)  9*  [  L.  ^  —j 

s  3T  (1  +  K)  (1  +  K) 

_  X 
aT  “  4ro 

v 

where 

\  =  wavelength  (meters) 

av  =  rms  clutter  velocity  (m/sec) 

So,  L  ,  the  loss  due  to  staggering  is  given  by 
s 


L  -  -  q  -  K)2  ,h_  16 _  [  2 

222  2  4*^ 

12T  tt  av  (1  +  K)  (1  +  K)  w 


(1  -  K)S 
12 


In  Figure  43  I  for  a  three-pulse  canceller  has  been  plotted  for  various 
stagger  ratios. 

Since  both  the  frequency  response  curve  and  the  KTI  improvement,  I, 
are  affected  by  varying  the  interpulse  period,  the  designs  developed  for  the 
unstaggered  case  are  not  necessarily  optimum  for  the  staggered  case  as  is 
illustrated  in  Figure  44.  The  results  shown  in  Figure  44  illustrate  the 
effects  of  prf  stagger  on  the  response  of  one  of  the  minimum-rms  error  filters 
developed  earlier;  in  fact  its  response  might  b.s  considered  less  desirable 
than  the  conventional  response  shown  in  Figure  42. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  Primary  Conclusions 

The  free-space  performance  of  the  EAR  system  may  be  seriously  affected 
by  changes  in  Doppler  frequency  of  the  target.  The  region  of  optimum  angle¬ 
tracking  accuracy  of  the  EAR  system  is  approximately  four  to  seven  kilometers, 
depending  upon  the  cross-section  and  physical  size  of  the  target.  Clutter 
cancellation  of  the  EAR  system  will  primarily  be  limited  by  an3log-to-digital 
conversion  quantization  noise. 

The  effect  on  system  performance  of  changes  in  Doppler  frequency  of  the 
target  may  be  minimized  by  using  samples  with  non-integral  sample  weights. 

The  design  procedures  set  forth  in  Section  V  produced  significantly  improved 
MTI  processors  for' the  unstaggered  prf  case,  but  very  little  has  been  Gone 
in  the  case  of  staggered  prf. 

System  perfor.  ance  may  be  seriously  limited  by  multipath  returns  when 
tracking  low-flying  targets.  These  tracking  errors  are  significant  for 
targets  flying  at  altitudes  less  than  approximately  1500  feet. 

Frequency  agility  may  be  of  appreciable  help  in  reducing  the  effect  of 
multipath  returns  on  system  performance.  The  effectiveness  of  frequency 
agility  in  reducing  multipath-induced  errors  is  a  strong  function  of  the 
frequency-agility  bandwidth  used,  the  target  heights,  and  the  antenna  height. 
The  400-MHz  bandwidth  which  may  be  possible  with  a  modified  EAR  would 
significantly  reduce  these  low-angle  tracking  errors,  particularly  for 
higher  antenna  heights. 

B.  Principal  Recommendations 

Based  on  the  information  studied  during  this  program,  the  following 
recommendations  are  presented  for  the  EAR  development  program. 

(1)  A  flexible  digital  MTI  filter  or  processor  should  be  incorporated 
into  the  system.  The  performance  obtainable  by  the  selection  of  sample  weights 
appears  necessary  if  acceptable  system  operation  over  a  wide  range  of  target 
Doppler  frequencies  is  to  be  achieved. 

(2)  Stagger  of  the  prf  should  be  incorporated  into  the  EAR  system  concept 
in  order  to  eliminate  "blind  speeds"  of  the  system,  and  investigation  of 
improved  MTI  filters  for  staggered  systems  should  be  initiated. 

(3)  Frequency  agility  should  be  incorporated  into  the  EAR  system,  both 
to  reduce  free-space  tracking  errors  due  to  target  glint,  and  to  reduce  low- 
angle  tracking  errors  due  to  multipath  signals.  Weighing  of  the  received 
signals  on  the  basis  of  the  sum  signal  amplitude  appears  desirable. 

Experimental  evaluation  of  the  EAR  shou/J  be  paralleled  by  careful 
theoretical  studies  as  an  aid  in  planning  significant  experiments  and  to 
assist  in  interpreting  results  of  those  experiments. 
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Appendix  A 


CONVENTIONAL  THREE-PULSE  MTI  FILTERS 


Consider  the  two-delay  filter  shown  in  Figure  A-l.  The  impulse  response 
of  this  configuration  is 

h(t)  =  6 (t)  -  2  6(t  -  T)  +  6 (t  -  2T) 

In  the  frequency  domain 

H(oj)  =  1  -  2eJ“c  +  e2** 

and 

G((u)  =  H(m)  •  H f  (u))  =  16  sin4 

It  is  interesting  to  examine  the  improvement  in  target-to-clutter  ratio 
when  this  canceller  is  used  as  an  MTI  radar  processor.  This  calculation  wi?.l 
be  made  assuming  a  Gaussian  frequency  distribution  of  clutter  power.  The 
total  received  clutter  power  will  be  denoted  by  C. 

Let  target  power  be  distributed  uniformly  in  frequency,  with  a  power 
of  T^.  The  input  target-to-clutter  ratio,  r^  is  then 

T. 

l 


The  output  signal  and  clutter  powers,  T  and  C  ,  must  be  computed,  and 
the  output  target-to-clutter  ratio,  rQ,  will  then  be 

T 

-  _2. 

*  o  C 

o 

The  improvement,  I,  if  any,  is  then 


Now, 


r  T  /C 
o  _  o  o 

r.  "  T./C 

i  i. 


f  n 

r 


(Clutter  power  spectral  density)  *  G(uO  df 


Preceding  page  blank 


89 


and 


If  is  often  customary  to  deiiue  cr  =  a  X/2,  where  X  =  wavelength  (meters). 

v  c 

Representative  values  of  c  are  given  in  Table  I,  and  Figure  A-2  shows 
I  plotted  as  a  function  of  a  ft*  values  of  f  =5  kHz  and  F  =  5.5  GHz,  where 
F  =  radar  frequency.  v  r 

Some  sample  calculations  of  the  improvement  obtainable  can  be  made  with 
the  use  of  Table  II. 

If  one  assumes 

O =  0.22  (wooded  hills,  20  kts) 

X  =  5.5  cm 

s 

f  =  5  x  103  Hz. 
then ,  r 

I  =  82.9  dB. 
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TABLE  I 


Source  of  Clutter 

Sparse  Woods 

Wooded  Hills 

Wooded  Hills 

Wooded  Hills 

Wooded  Hills 

Sea  Echo 

Sea  Echo 

Sea  Echo 

Sea  Echo 

Chaff 

Chaff 

Chaff 

3a in  Clouds 
Rain  Clouds 


SUMMARY  OF  STANDARD  DEVIATIONS  OF  THE  CLUTTER  SPECTRUM  [52] 

L Cjutter  Wind  Speed  (knots)  av  (m/sec) 


Calm 

0.017 

10 

0.04 

20 

0.22 

25 

0.12 

40 

0.32 

— 

0.7 

-- 

0.75  -  1.0 

8-20 

0.46  -  1.1 

Windy 

0.89 

— 

0.37  -  0.91 

25 

1.2 

-- 

1.1 

-- 

1.8  -4.0 

■k-ruu' 


Appendix  B 


LIMITATION  IN  IMPROVEMENT  FACTOR  CAUSED  BY 
ANALOG-TO-DIGITAL  CONVERSION 


Since  radar  systems  using  digital  processing  must  operate  with  a  finite 
number  of  bits,  errors  are  introduced  by  the  process  of  analog-to-digital 
(A/D)  quantization;  these  errors  impose  a  limit  c.n  the  Mil  improvement  that 
can  be  obtained.  In  a  system  for  which  E  is  the  value  of  the  least-signifi¬ 
cant  bit,  there  is  a  possible  error  +E  / 2° in  each  conversion.  This  error 
has  a  uniform  probability  density  if  tRe  least-significant  bit  is  uncorrelated. 

Let  the  error  in  quantization  be  e,  and 
e  =  0 


0E/2  0 

a =  f  °  (x  /E  )dx  =  E/12, 

6  -E  /2 
o 

2 

where  is  the  quantization  noise  power. 

In  a  nonrecursive  filter  such  as  those  studied  earlier,  the  noise  power 
out  of  the  filter  due  to  quantization,  P  ,  is 


E 

p  S  S-. 

n  12 


Vi 

Z  ■ 


It  has  been  shown  [53]  that  if  the  filter  output  is  from  -1  to  +1,  the 
pulse-to-pulse  signal  deviation  due  to  quantizer  noise  is 


(2‘*  -  1  )  /O' 

where  M  is  the  number  of  A/D  converter  bits. 


The  limit  this  imposes  on  the  improvement  factor,  I,  is 
I  *  20  log  (2  -  1)  A. 5 

Since  in  a  two-channel  Mil  there  are  two  independent  errors,  then  at 
the  output, 

I  =  20  log  (2h  -I)  ATT  . 

For  a  nine-bit  A/D  converter,  this  is  52.9  dB,  and  for  a  nine-bit  plus 
sign  conversion  (ten  bits)  this  is  58.9  dB. 
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RANGE  TRACKING 


In  conventional  radar  systems  using  analog  signal  processing,  the 
so-called  split-gate  range  tracker  is  often  used  to  track  a  desired  target 
in  range.  In  a  radar  using  digital  signal  processing,  however,  the  conventional 
split-gate  range  tracker  cannot  be  used,  since  received  signals  are  sampled 
at  discrete  points.  Nevertheless,  the  information  obtained  from  these  discrete 
samples  may  be  used  to  perform  a  range -tracking  function. 

Consider  a  representative  received  pulse  (figure  C-l)  which  is  sampled 
at  points  spaced  A  units  of  time  apart.  The  current  estimate  of  pulse 
position  we  designate  by  t  and  the  pulse  amplitude  at  t  by  E  .  Associated 
with  E  are  two  other  samples,  E  and  E^,  which  occur  at  times  t  -A  and  t  +A 
respectively.  These  three  samples  may  now  be  used  to  calculate  the  error°in 
range,  e  ,  between  t  and  some  point  on  the  received  pulse.  A  number  of 
computational  algorithms  are  possible;  one  which  appears  very  similar  to  the 
analog  split-gate  range-tracking  scheme  is 


where  K  is  a  porportionality  constant.  Another  algorithm  which  has  been 
proposed  by  the  General  Electric  Company  is 


£*+£o 


A  brief  analysis  of  the  performance  of  these  two  range  tracking  algorithms 
has  been  undertaken.  There  are  a  number  of  variables  in  such  an  analysis 
that  affect  system  performance,  for  example:  the  particular  algorithm,  the 
transmitted  pulse  shape,  the  receiver  bandpass,  and  the  sample  spacing.  In 
order  to  simplify  the  analysis,  a  rectangular  transmitted  pulse  and  a  Gaussian 
filter  response  were  assumed. 

The  first  step  in  the  analysis  was  to  calculate  the  received  pulse  shape 
for  a  rectangular  transmitted  pulse  and  the  appropriate  receiver  bandpass. 

A  Gaussian  filter  was  assumed  since  this  closely  approximates  the  response 
which  will  probably  be  implemented  in  the  EAR  system. 


The  frequency  response  of  a  Gaussian  filter  is  given  by  [54] 


H(tii)  -  Koe-°-^\y 
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where , 


H(to)  “  voltage  output  of  the  filter, 

A  =  a  scale  factor, 
ui  =  angular  frequency  (2rrf) 

u>c  =  3-dB  video  bandwidth  of  the  filter  (in  radians) , 

and  the  filter  delay  has  arbitrarily  been  set  equal  to  zero.  In  the  above 
expression,  u>  >  the  3-dB  video  bandwidth,  is  equivalent  to  half  the  3-dB 
bandwidth  at  £he  IF  frequency. 


The  response  of  such  a  Gaussian  filter,  a(t),  to  a  unit  step  function 
is  [55] 


n  /  0)  t  \ 

-  f  i1  +  er£  aksr)  • 


where 


9  x  2 

erf  x  =  Jo  e’y  dy‘ 


Therefore,  the  filter  response,  f(t),  to  a  rectangular  pulse  of  width,  T, 
is  given  by 

A  ,  a)  t  .  A  ,  u>  (t  -  t) 

£(t)  -  f  (l  +  «£  5^)  -  f  (l  +  )  , 


or 


f(t)  =  ~  |erf  (0.845  u^t)  -  erf  [0 . 845 (t  -  t)  u)c]} 


This  information  was  used  to  calculate  sample  values  £  ,  £  ,  and  £  . 

The  range  error  was  then  calculated  using  each  algorithm.  lesuJts  are  plotted 
as  functions  of  e  /K  versus  normalized  time  error  t  /t.  The  condition 
t  h  -  0  corresponds  to  the  on-target  sample  coinciSing  with  the  leading 
e9ge  of  the  rectangular  pulse  input  to  the  filter,  while  t  /t  -  1  corresponds 
to  the  trailing  edge.  Figure  C-2  shows  e  /K  as  a  function  of  t  h  for  the  case 
^3dBT  ”  *  ^3dB  *s  t*ie  ^  3-dB  bandwidth)  Ind  sample  spacing  A/t°=  1.  Figure 
C-3  is  the  same  type  of  presentation  for  f~  -  2.  The  actual  performance 
of  the  EAR  system  probably  lies  somewhere  Between  these  two  cases. 

The  performance  of  the  two  tracking  algorithms  in  the  presence  of  receiver 
noise  is  an  important  system  design  consideration  While  an  exact  analysis 
is  beyond  the  scope  of  this  project,  an  approximate  analysis  gives  useful. 


insight  into  the  performance  of  the  two  algorithms. 


and 


Assume  n^,  n  ,  and  are  the  noise  voltages  associated  with  E  , 
?£,  respectively.  Then  for  the  conventional  tracking  algorithm6 


2  , 


e 


r 


=  K 


L- 


(Ei  +  n.)  -  (S  +  n  ) 


E  +  n 
o  o 


Assume  the  signal-to-noise  ratio  is  greater  than  one,  such  that 
E  +  n  »  E 

O  O  0 

Then 


e  « 
r 


K 


E 

e 


+  n*‘ 


E 

o 


n 


e 


If  n£  and  n  are  uncorrelated,  and  their  noise  powers  are  both  equal  to  n, 
and  che  range  tracker  is  on  target  (E^  -  E  =  0) ,  then  the  rms  range  error 
due  to  noise,  is  given  by 


If  S/N  is  the  signal-to-noise  power  ratio,  then 


a 


rn 


=  K 


[7T  1 

j/s7n_ 


While  strictly  true  only  for  large  S/N,  the  general  prediction,  that 
range  errors  are  directly  related  to  K  and  inversely  related  to  the  signal- 
to-noise  ratio,  appears  to  be  indicative  of  the  performance  for  other  cases 
as  well. 


Figures  C-2  and  C-3  show  that  e/K  is  a  much  more  nonlinear  function  of 
range  error  for  the  conventional  algorithm  than  for  the  General  Electric 
algorithm.  They  also  show  that  the  conventional  algorithm  is  relatively 
insensitive  to  small  variations  in  position  about  the  true  pulse  location, 
particularly  for  the  case  f  t  =  2,  as  evidenced  by  the  small  slope  of 
e  / K  at  t  /t  =  0.5.  Because  of  the  nonlinearity  of  the  conventional  algorithm, 
iE  would  appear  that  any  assumed  value  of  K  for  tMs  system  would  be  only 
approximately  correct  over  the  region  0  £  t  /t  £  1;  it  appears  that  this 
average  K  would  closely  approximate  the  value  of  K  for  the  General  Electric 
algorithm,  so  the  noise  performance  would  be  approximately  the  same  for  the 
two  cases.  These  results  indicate  substantial  advantages  in  using  the  algorithm 
proposed  by  General  Electric.  Realistic  values  of  K  for  the  two  algorithms 
are  comparable,  indicating  comparable  noise  performance.  The  greater  linearit} 
of  the  General  Electric  algorithm  and  the  fact  that  it  remains  bounded  for  large 
errors  offer  substantial  advantages.  Partially  offsetting  this  advantage  is 
the  fact  that  the  conventional  algorithm  tracks  with  a  sample  on  the  center  of 
the  pulse,  thus  providing  the  maximum  signal-to-noise  ratio. 
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Appendix  D 


EFFECTS  OF  FREQUENCY  AGILITY  ON  BLIND  SPEEDS 


Changing  the  carrier  frequency  of  MTI  radar  has  no  effect  on  the 
frequency  response  of  the  filters,  but  the  target  blind  speeds,  V  ,  are 
shifted  slightly,  since  ° 

VB  =  N  ~‘2^P"£  (knots)  N  =  1,  2,  .  .  .  . 
prf  =  pulse  repetition  frequency  (Hz) 

F  =  radar  frequency  (GHz) . 

In  a  frequency-agile  radar  the  blind  speeds  vary  from  pulse  to  pulse, 
yielding  a  velocity  response  that  is  the  sum  of  the  output  at  each  frequency. 

Figures  D-l  and  D-2  show  the  response  for  a  frequency-agile  radar 
operating  with  two  frequency  step  sizes.  Figure  D-l  is  for  50-MHz  step  size 
and  Figure  D-2  is  for  31. 25 -MHz  step  size.  The  starting  frequency  is  5.3  GHz, 
and  the  stopping  frequency  is  5.6  GHz  for  Figure  D-l  and  5.8  GHz  for  Figure  D-2. 

It  may  be  noted  that  a  marked  change  is  made  in  the  velocity  response 
without  benefit  of  pulse  staggering,  but  approximately  20  dB  null.,  are  present 
at  the  lower  blind  speeds. 
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Appendix  E 


COMPUTER  PROGRAM 


In  this  appendix  is  presented  the  FORTRAN  V  computer  program  used  in  the 
computation  of  multipath  errors  in  the  EAR  system.  The  subroutines  listed 
are  representative  of  those  employed  in  the  generating  data  for  this  report; 
however,  various  slightly  modified  versions  were  used  as  necessary  to  produce 
special  outputs. 
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' _ 91  meN5i O.N_I BU^ (500 QbVTl TL£ll2),_HT I TLE U l) _ 

call  PLors(nurd), 5ooo*2) 

_ WRIT-(6»90q) _  _  _ _ 

900  FORMAT <  *QNO  PREDICTION  FILTER*  FREQUENCY  AGILE*) 

s _ OlME^SION  NARRAY(5000) .X(IOOO) .ERRFT(inOO) _ 

DIMENSION  SLOGG(IOOO) 

! _ PIZATAN2  ( 0*i-l  .J _ 

RO=l8o,/P'l 

lOo  FOR’lAj ( ) 

_ 101  FORMAT ( » OEnTER  HA»HT*RSTART*V*pRF*FrEQ(GHZ) *RHO* ) _ 

102  FORM AT {*  3yE ' ) 

_ Pl02=PI/2 . _ 

IOSTRTz-200 

_ I9.S.T0.P..=2QQ _ 

ILAST=900 

_ ICNT=4Q0 _ 

THSCANsO. 

_ WRITEtG^BOQ)  _ 

800  F0RmAT(*6ENTER  THANT»OlEVEL*) 

_ REA3<5* 100)  THANTtOLEV _ 

RHOzO. 

_ AKS=. 85934 _ 

10  WRIT- (6* 101 ) ” 

_ read(5*ioo.endz999*errz999)  ha,ht*R5Tart *v*prf*freq»rhq 

WRI T- ( 6*  94o ) 

940  FORMAT (*0£nT£R  fSTART*FSTQP,FStEP*) _ 

READ <5* 100)  FSTART,FSTOP»FSTEP 

_ RANGEz_RSTA&T*6000, _ 

THTS=ASIM((HT-HA) /RANGE) 

_ IHSCAn=-( THANT-THT5*R0) _ 

DELTATZ450/PRF 

_ I=Qj _ 

IFRrOziNTt (FST0P-F5TART)/FSTEP)+1 

9.18 _ CONTlNgE _ 

TOzT 

_ FREO=FSTARX-crSTE=' _ 

.  00  20  lzi,3' 

_ FHQD-MQ0 ( 45O  > IFREQ ) _ 

FRE0=FREQ*F'10)*rSTEP 
ALAM- , 9843/pRE0 
thscanzthscan+anger'r 

CALL  ERRO ( HA, THANT* HT. RANGE *ALAM*A!tpERR*V*RHO*THT»THSCAN* 
C  FSTART,FSTOP,FSTEP*PrF*OLEV*9LOGG(I)) 

_ T=TfO£LTAT _ _ 

XXZTHSCAN4-THANT 

_ X(I)=XX+ANSERR _ _ 

ERRzXx-THT  ‘ 

•  _ ERRFT(I)zRaNGE*ERR/RO 


T*fREO* 


P I.311  re  K-l.  Ma i r.  Computer  program. 
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Y=R AW5£*C0S ( THT /RO ) -V*6 . *6t L T A  f/SVb 
_ RAno£=SQRT(Y*«2+HT-HA)**2) 

20  cOn't  In  i  j  £  _____ - 

_ LCNT-200 

LCf =LCNT-3 
_ 00_4C  I=1»lCT 

FRE Q-FRE G+F MOO *F STEP  “  - - 

ALAM=.9943/-R£Q 
THSCA,NJ:TH5C  AN*ASi  j£o.R 

_ CALL_ERRO_(HA,THANT.HT»RANSE#  ALAVI,  ANr,ERR,V»r>HO.THT#THSCA*nT#FREO» 

c  F5TART»FST0P»FSTE?YPrF#DLEV»5L0GG( 1+3) j ” 

_ IrI+°ELTAT 

XX=fHANT<-THSCAN  “ 

_ Xj I +3) =XXf ANGERR 

ERR=XX-THT  — — — 

_ ERRFT ( I f 3) =RANGE*ERR/RQ 

Y=RA'NGE*COS(THT/RO}-V*6V*OELfAT/3,6 


RANGE=S0RT(Y**2+(HT-HA)**2) 

_ FREf)l=PREQftOOOOOl  _ 

iF(FREOl.GE.FSTOP)  pRSQ=FSTART-FSTEP 

...40. —.CONTINUE  _ _  _ _ _ _ _ 

CALL  DIST(NARRAY.ERRFT,IOSTRT*IDSTOP.LCNT»ICNT,IOVER* HINDER) 

_ WRIiLCfe 1 9231 _ 

929  FORMAT! 'OEnTER  1  TO  PLOT  DATA*) 

_ S- A D  <  5  r  1 0 0 )  NPLOT  _ _ _ _ 

NCNT=NCNT*1 

_ IFlNpLOT.N£.l)  SO  TO  923  . _ _ 

IF(mCnt.GT.I)  GO  TO  929 

_ DO  777  1  =  1. IP _ 

VTITLE(I)=*  • 

777 _ CONTINUE  __  _  _ 

VTI Tl»S ( 1 )  =  •  ERROR  ♦*  '  * . " 

_ ymLE.(2)=*("EET5'_ _ _ _  _ _ 

VL3nD=-400. 

_ yua^DsMQ.. _ 

926  FORMAT! 12A6) 

__  WRIT4(6f927)  __  .  ..  _ _ _ _ 

927  FORMAT ( *  OENTER  TITLE  FOR  HOR  AXIS') 

.  REA}(5*926)  HTI TLE  _  _  _  _  _ 

929  HL3'|D  =  T0 

_ HU3'P=J _ 

CALL  SI MPL (ERRFT »LCNT*vL3ND»VU3NDHL9nd»HURMD»VTITlE*HTITLE) 

923  CONTlNJE  _ _  _  _  _ 

MM=MV'  +  LCNT 

_ WRIT-(6»920)  RANGE, vm 

920  FORMAT! 'ORANGE  =,,P7,1,•  FEST»*»I5,'  ITERATIONS') 

_ WRIJE.!  6.t  9 1 9J _ 

919  FORMAT! 'OENTER  1  TO  CONTINUE  *) 

_ RE AO <S, 100)  NCONT  _  _ 

IFINCONT.EO. 1 )  SO  TO  919 

999  CONTINUE  _  _  __ 

CALL  PLOT(10,»0.'909) 

_ EM _ 


F  1  "lire  r  - 1 . 


i  font  *i  •  ) 


COMPLEX  PUmC TION  GAIN ( tH , SCAMTH>N) _ 

THSC7'N=AiS(SCANfH) 

complex  t»arg _ 

Pr02=ATAM2(l,  *0,1 

PI=2+PI02  _  _ 

'  RD=lBd,/Pf'"' 

_TH2£«0=.73_  +  « 27* ( THSCAM/40« ) _ 

F=1,5T  •  5*  ( THSCA,N/40* ) 
_THM=(TH-rH7-PnwRD 
THP=(TH+THZER0)/R3 

UM=-*33.9*PI*SIM(THM)  _  _ 

UP=?*33,9*PI*SIN(THP) 

_TM£5 *PJ 02*72 . 8* COS ( UM ) / (PI02»»2-UM**2) 
TP=.-*PI02*72ra*C0S(UP)/(PI02**2-UPt*2) 

ARG=CMPLX(0. »PI02) _  _ 

IF(N.EO.'l)  T=TM+TP 

IF  (N.£0._2L  T=(TM-TP)»CEXP(ARG) _ 

GAI<|=Y 

RETURN) _ 

END 


Figure  E-2.  Subroutine  GAIN 
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_ _ SUSROut  I  NE_  0 1  ST  <  NARR  A  Y ,  Xj  I  START?  I SJQP  tj  COUNT?  ICNJ  tI  qVER?IUNQER) 

DIMENSION  NARRAY(5nOO) 

IQv-sRso _ _ 

iund£r=0 

_ OIM£NSION_x1U _ 

K=0 

iAQ-fOP'^^ODISTRIBUTlON  RANGES  £c?OM  T0«?l6t*.*) _ 

“101  FORMAT ( *o there  WERE  r?  1 4?  > 'points  ThAT  v:ERE  out  OF  range  of  the  di 

_ oSTRIBUTION  array?) _ 

1 02  FORwA.t ( ?  o  MEOIAN  IS* » 25) 

JOa„£PR«AT{  !0  ?.,!<♦?  •  POINTS..  WERE  UNDER. POINTS  WERE  OVER.*) _ 

M3 I NS= I STOP- I ST ART 

_ if  (  I  stop,  l=;,istartj_s.o  _to  99? _ 

IF (NBInS, 3r, 5000 )  GO  TO  999 

_ IMl£15Iflg _ _ _ 

irOP=ISTARj 

_ PiLlQ.  L~1*N?I/4S _ 

lo  nArrAY(I)=o 

_ POj>0_  I  =  1,1  COUNT. _ 

IX=lNT{X(I)+0.5)-ISTART 

IF(IX.GT.ITOP)  itop=ix 

_ IF(IX,LV.,I30I)  JROTSIX  _ _ 

IF(IX,ST,.N3INS)  I0vER-I0vER+3 

_ !F.lI>',LT.t.l)„XUNDER=IUNDER+l _ 

IF((IX.GT.N3INS),OR.«2X.LT.l))  GO  TO  30 


|  GO  JO  20 

|_3q  K=X  +  I _ 

1  2q  CONTINUE 


icnt=icount-k 

_ ,I30T=I30T  +  ISIART _ 

lfoo=ITOP*ISTART 

_ WRIT?  <  6  ?i  0  0  )„  I  30.Tp.II0S _ 

SUM;0, 

_  QO  UP  ,!  =  I?N3INS _ 

K=I 

_ SUMrNARRAY (I ) /FLOAT ( ICNT) +5UM 

IF<sOM,3E.«5)  GO  TO  999 
4q  continue _ _ 

99g  CONTINUE 

_ X  =<  » I  ST  ART _ _ 

RETURN 

END _ _ 


*>'■.  Jubr-’i".  i no 
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__  SUBROUTINE  STNQEV(X»ICOUNT)  _  _ 

J INSMS I ON  X<1) 

_ S=0.t _ 

SQ=0. 

_ 00  10  i=t.»icogRT. _ . _ 

S=S+X(I) 

xsa=x(i)**2 _ _ _ 

fi'J=S0+XS3 

,10  continue  _ _ „ 

XAV3=S/IC0uNT 

.  .  XSQ/\Vr,=S0/lC0UNT  ..  _ _ 

S0-S0RT(XS0AVG-XAV6**2) 

WRIT£(6*100>  XAVfi.  . . . 

100  FORMAT!  »0  AVERAGE  ISSF6.1) 

_ WR I  TE(6»20  o  1_SQ _ 

200  FORMAT! *0  STANDARD  DEVIATION  IS’'F6,1> 

END  ' 


Figure  E-4 .  Subroutine  STNDEV 
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SUBROUTINE  ER30<HA>THAMT,HT«RfeNGEtALAMtANSERR»V>RHQtTHTtTHSCAN,T» 
C  ~  ~  FREQ»FsTART»cS‘*OPiFsTEP»PRF*OLEV»SLOG1)  '  r 

_ COvpLE_x_TA,Ta.T5*T«»T9.T10#Tll,Tl2 

DIMENSION  XAMStfOO)  OL03UOO)  »SL03(jOO) 

_ COMPLEX  GAl*J;EXP»RATIO 

COMPLEX  T'i,T8»T6»f7»S*b 

„ _ AKS=.  845*.  335*  iTMSC  AN/30.) 

OELfATsi/PRF 

PI=ATAn2(0.>«1.) _ 

ro=i6o./pi 

_ X=0. 

I F TF Peg • G E , FST 0 P -FST EP )  FREOsFsTART-FStEP 
_  DO  10  I=1#4P 

IMOD=MOOa,3)  *  ~  " 

_  IF(I^OD.EO.O)  FREQs^REO+FSTEP 

IF (7^00. £3, 0)  ICNT=1 

_ ifjimod.ne.o)  icntsicnt+i _ 

ALAMP, 9843/FREQ . ~ 

THT= A  SIN!  <  ViT-HA) /RANGE) 

_ THT=Th,T*R3  ___ _ _ 

YYssGrt <R AnSE**8-CHT»HA > **2) 

_ THIM=-ATAN2(  (HA+HTJ.,  Yv  }  _ _ 

THIM=TH1M*ro 

_ Y=YYrVOELTAT/,6 _ _ 

RAN3L=SQRT<Y**2+(HT-HAV**2> 

_ _ 0ELTAr=4,*pUhA*(HA+HT,'/(ALAW*rANGE)_  _ _ 

TH3EAM=THAnT+THSCAM 

_ TR1=Tht-H3EAM _ _ _ _ _ _ 

Tri2=THIM-TH3EAM 

EXP=CMPLX(n.>->OELTAR)  _ 

TI=SAIN(TH1#THSCAN,1> 

_ ia=GAI'j(Td2*THSCAN,n  _ _ 

T6=5AINITH1»THSCAN,25 

_ T7=5AlN(TH2fTHSCAN,2}  _  _ _ _ 

T8=T64RCS(TH1) 

_ L9=ilH0*c  E  X  RiEX£l±I<?  iR.CS  LlltU _ 

T10=RHO*CExP(-XP)*T7*RC.SfTd2) 

_ Tll=Pd0*RH0*CEXPU*EXDWT7^RCSf,TH2)  _ _ _ _ 

D=Tfl+T9  +  Uo  +  m 

_  DMA3=CA3SO)  _  _ _ _ 

IFCjVa3.LT. IE-10)  ?VAG=1E-10 

_ ?l.OO±D_=2  0  *  ALPGJLmvJi5J _ 

T3=Tl*RCS(rU> 

T4=RHo*cEXP(EXP) *T1*RCS« Ml)  _  _  _ 

T5=R^0*CEXP(EXP)*T2*RC9(T'-I2) 

T12=RH0*RH0*CEXP(2*EXP)*T2oP!'  WTH2)  __ _ 

S=T3+T4+T5+T12 


F i *;»u >"  K-').  ->ul>rout.ine  KRRO. 
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FUNCTION  RC5(TH) 
OATA  30/57.29578/ 
THRsTH/RO 

„TH  R  X- 1 00  *TH/PO. _ 

mC0S{THR)**2 

V2=C0S(THRi) 

\“J-*  ♦ 

RETURN _ 

END 


Figure  E-6 


Subroutine  RCS. 


APPENDIX  F 


FREQUENCY  RESPONSE  OF  NONRECURSIVE  DIGITAL  FILTERS 


The  nonrecursive  digital  filter  is  illustrated  in  figure  F-l.  This 
processor  is  typified  by  simplicity  ot  analysis,  economical  hardware  imple¬ 
mentation,  and  short  settling  time. 

The  output  of  the  filter  is  formed  from  a  weighted  sum  of  the  previous 
N  +  1  input  samples,  where  N  is  the  number  of  delay  elements.  It  follows 
directly  from  Figure  F-l  that  h(t),  the  impulse  response  of  the  filter,  is 
given  by 

h(t)  =  A  5(t)  +  k  6(t-T)  +  A2  <5(t-2T)  +  ...  +  A..  (t-NT) 

N  a) 

-  k=0 

where  6(t)  represents  a  unit  impulse,  and  Afi  is  the  weight  of  the  Nth  sample. 

To  compute  the  complex  frequency  response  of  the  network,  H(oj),  the 
Fourier  Transform  is  applied  to  (1),  yielding 

«<“>  ■  R=0  \  e'jk"T  '  (2) 

The  power  response  of  the  filter,  G(u),  is  then 

N 

G (to)  =  H(u>)  .  H  ( u> )  =  C  +2  C  cos(qwT),  (3) 

o  q=l  q 

N-q 

£ 

where  C  =  „  A.  A.  .  0  <  q  <  N. 

q  3=0  k  k+q  —  ^  — 

Equation  3  will  be  recognized  as  a  truncated  Fourier  Series,  whose 
coefficients  are  simple  combinations  of  the  filter  coefficients. 
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